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PREFACE 


The  model  study  of  the  velocities  induced  by  commercial  navigation  was 
authorized  by  the  US  Army  Engineer  Division,  Ohio  River  (ORD) ,  at  the  request 
of  the  US  Army  Engineer  District,  Louisville  (ORL) . 

The  study  was  conducted  during  the  period  July  1987  to  September  1988 
in  the  Hydraulics  Laboratory  of  the  US  Army  Engineer  Waterways  Experiment  Sta¬ 
tion  (WES)  under  the  direction  of  Messrs.  F.  A.  Herrmann,  Jr.,  Chief,  Hydrau¬ 
lics  Laboratory;  R.  A.  Sager,  Assistant  Chief,  Hydraulics  Laboratory;  and 
G.  A.  Pickering,  Chief,  Hydraulic  Structures  Division  (HSD) ,  Hydraulics  Labo¬ 
ratory.  The  tests  were  conducted  by  Dr.  S.  T.  Maynord,  Project  Engineer, 
Spillways  and  Channels  Branch  (SCB),  HSD,  and  Messrs.  D.  White  and  J.  Hilbun, 
SCB,  under  the  direct  supervision  of  Mr.  N.  R.  Oswalt,  Chief,  SCB.  The  report 
was  written  by  Dr.  Maynord  and  edited  by  Mrs.  M.  C.  Gay,  Information 
Technology  Laboratory,  WES. 

During  the  course  of  the  investigation,  Messrs.  L.  Richardson,  ORD,  and 
D.  Beatty,  T.  Siemsen,  J.  Kleckner,  J.  Baker,  and  B.  Vessels  of  ORL  visited 
WES  to  observe  tests  and  discuss  test  results. 

Commander  and  Director  of  WES  during  preparation  of  this  report  was 
COL  Larry  B.  Fulton,  EN.  Technical  Director  was  Dr.  Robert  W.  Whalin. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non— SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ £* _ 

0.01745329 

0.3048 


To  Obtain 
radians 
metres 


horsepower  (550  foot-  745.6999  watts 

pounds  (force)  per  second) 

inches  2.54  centimetres 


miles  (US  statute)  1.609347  kilometres 

square  feet  0.09290304  square  metres 
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VELOCITIES  INDUCED  BY  COMMERCIAL  NAVIGATION 

PART  I :  INTRODUCTION 

Background 

1.  As  a  vessel  moves  through  a  body  of  water,  it  creates  disturbances 
in  the  form  of  altered  water  levels  such  as  waves  or  drawdown  and  altered 
velocity  patterns  such  as  increased  velocities  and  turbulence  intensity. 

These  disturbances  vary  in  intensity  with  such  factors  as  vessel  type,  speed, 
size,  and  location;  and  channel  size,  shape,  and  ambient  conditions.  The 
intensity  varies  significantly  with  distance  from  the  vessel.  In  large  water¬ 
ways  this  disturbance  may  be  greatly  diminished  by  the  time  it  reaches  the 
boundaries  of  the  bed  and  bank.  The  ability  to  predict  the  intensity  of  these 
disturbances,  both  near  and  far  removed  from  the  vessel,  is  needed  to  better 
assess  the  physical  impacts  of  navigation  on  environmental,  channel  stability, 
and  other  concerns  within  a  waterway.  At  present,  tools  to  predict  the  inten¬ 
sity  of  these  disturbances  are  quite  limited.  The  extensive  body  of  European 
literature,  some  of  which  will  be  referenced  herein,  is  primarily  applicable 
to  channels  having  a  relatively  small  blockage  ratio  defined  as 

^  waterway  cross-sectional  area* _ 

submerged  cross-section  area  of  vessel  ^  ’ 

Navigation  in  many  US  waterways  exists  on  large  rivers  that  have  large  block¬ 
age  ratios.  Information  is  needed  on  navigation  effects  in  waterways  having  a 
wide  range  of  blockage  ratios. 

Purpose  of  Study  and  Outline 

2.  The  purpose  of  this  study  was  to  improve  the  understanding  of  flow 
patterns  induced  by  moving  tows  and  to  develop  tools  for  predicting  the  in¬ 
tensity  of  disturbance  for  navigation  and  waterway  sizes  typical  of  those 


*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  defined 
in  the  Notation  (Appendix  B) . 
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found  in  the  United  States,  particularly  in  the  Ohio  River  basin.  This  study 
focused  on  the  velocities  induced  by  commercial  shallow-draft  navigation. 
Velocities  resulting  from  both  the  propeller  jet  and  from  the  displacement 
effects  of  the  vessel  were  evaluated  in  this  study. 

3.  The  study  will  be  documented  in  the  following  manner: 

a.  Literature  search. 

b.  Flow  visualization  with  a  1: 20-scale  model. 

c.  Flow  visualization  and  surface  current  patterns  with  the  Olmsted 
Locks  and  Dam  model. 

d.  1: 20-scale  physical  model  study  of  bottom  velocities  near  the 
path  of  a  moving  tow. 

e.  Prototype  data  collected  by  the  US  Army  Engineer  District 
(USAED) ,  Louisville. 

f.  Development  of  predictive  relations. 

Unless  stated  otherwise,  all  quantities  are  given  in  prototype  values. 


5 


PART  II:  LITERATURE  SEARCH 


4.  The  literature  search  will  be  broken  down  into  the  following  areas: 

a.  Propeller  jet  studies. 

b.  Displacement  and  return  velocity  studies. 

c.  Field  measurements  of  vessel— induced  velocity. 

Propeller  Jet  Velocity  Studies 

5.  For  the  case  of  maneuvering  navigation  (vessel  speed  V  -  0) ,  the 
distribution  of  velocities  has  been  studied  by  several  investigators  including 
Fuehrer  and  Romisch  (1977),  Blaauw  and  van  de  Kaa  (1978),  Bergh  (1981),  Louis 
Berger  and  Associates  (1981),  and  Prosser  (1986).  The  assumption  of  maneuver¬ 
ing  navigation  greatly  simplifies  the  problem  because  wake  effects  are  elimi¬ 
nated.  Many  of  these  studies  involved  prediction  of  bottom  velocities  for  use 
in  determining  riprap  size  or  scour  depth.  For  vessels  underway  (V  *  0) , 
Fuehrer,  Romisch,  and  Engelke  (1981)  stated,  "Consequently  a  marked  reduction 
of  bottom  velocity  occurs.  Furthermore  the  maximum  bottom  velocity  takes 
place  In  an  ever-increasing  distance  behind  the  ship."  Gucinski  (1982)  con¬ 
ducted  laboratory  and  field  tests  and  determined  that  the  moving  vessel  resus¬ 
pends  less  sediment  than  the  stationary  vessel  having  the  same  propeller 
speed.  Velocity  measurements  demonstrated  that  the  propeller  jet  extends  to  a 
greater  depth  with  the  stationary  vessel.  Gucinski 's  results  also  showed  that 
turbulence  intensity  at  the  bed  was  decreased  with  the  higher  tow  speeds. 
Schale  (1977)  found  that  "the  propeller  jet  of  moving  freight  motorships,  even 
with  a  high  propeller  loading,  never  comes  in  contact  with  the  canal  bottom 
but  always  rises  along  the  shortest  path  to  the  surface  of  the  water."  Schale 
states  that  the  propeller  jet  strikes  the  channel  bottom  under  only  the 
following  conditions: 

a.  Startup  from  a  stationary  condition. 

b.  Whenever  the  water  depth/draft  ratio  is  less  than  1.2. 

c.  Maneuvering  with  hard  rudder. 

Schale 's  observation  of  the  jet  rising  to  the  surface  is  consistent  with  the 
findings  of  Maxwell  and  Pazwash  (1973)  for  shallow,  submerged,  axisymmetric 
jets.  Schale  (1977)  also  states  that  "the  propeller  jet  itself  is  undetect¬ 
able  in  the  measurement  plane  lying  0.5  m  above  the  bottom  while  in  fact  there 
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even  prevails  there  a  counterflow  in  the  propeller  direction.  Thus  the  at¬ 
tacks  upon  the  bottom  material  are  attributable  only  to  forces  arising  from 
displacement  flow  and  from  the  therewith-integrated  energy  field  of  the  ship's 
drive."  This  "energy  field  of  the  ship's  drive"  is  believed  to  be  the  flow 
entering  the  propellers.  Prosser  (1986)  presents  a  plot  of  bottom  velocity 
due  to  propeller  inflow  versus  blade  tip  clearance  above  bed. 

6.  Before  discussing  the  applicable  propeller  jet  studies,  the  basic 
relations  used  in  many  propeller  jet  studies  will  be  developed.  Propeller 
jets  are  often  considered  to  be  similar  to  submerged  jets  such  that  the  same 
equations  can  be  used.  Albertson  et  al.  (1950)  defined  the  velocity  distribu¬ 
tion  for  a  submerged  jet  discharging  into  an  infinite  fluid.  The  following 
assumptions  apply: 

a.  Pressure  is  hydrostatic  throughout  flow. 

b.  Diffusion  is  dynamically  similar  under  all  conditions. 

c.  Longitudinal  component  of  velocity  varies  according  to  normal  or 
Gaussian  probability  function. 

Unless  otherwise  stated,  equations  are  presented  in  a  form  applicable  to  any 
set  of  units.  The  equations  for  flow  through  an  orifice  are  as  follows: 


a.  For  the  zone  of  flow  establishment: 


y  =  exP 


-  _  ' 

r  +  Cx  - 

2' 

2(Cx)z 

(2) 


where 

Vx  =  velocity  in  x-direction  at  coordinates  x,r 
V0  =  orifice  velocity  at  outlet 
r  =  radial  distance  from  center  of  outlet 
C  =  coefficient 

x  =  distance  from  outlet  measured  along  jet  axis 
D0  =  orifice  diameter 
b.  For  the  zone  of  established  flow: 

v(x)m«  1  Dc 

=  7C  IT 
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where  V(x)max  is  the  velocity  in  x-direction  at  x,r  -  0 
and 


Vx 

v^y 


max 


exp 


(4) 


Based  on  the  experimental  results  of  Albertson  et  al .  (1950),  the  flow  becomes 
established  at  x/Dc  —  6.2  and  C  =  0-.  081  for  the  orifice  discharging  into 
an  infinite  fluid.  Equation  3  becomes 


V(x) 


max 


*  .6.2 


(5) 


It  is  important  to  note  that  propeller  jets  behind  moving  vessels  differ  from 
the  conditions  addressed  by  Albertson  et  al .  (1950)  in  the  following  ways: 

a.  The  channel  bottom  and  water  surface  inhibit  jet  spreading. 

b.  A  moving  jet  is  discharging  into  a  moving  flow  field. 

c.  The  propeller  jet  has  a  radial  component  of  velocity. 

d.  The  rudder  splits  the  jet  into  two  jets. 

e.  The  Kort  nozzle  and  open  wheel  are  different  from  an  orifice. 
The  following  paragraphs  summarize  the  five  methods  found  in  the  literature 
for  propeller  jet  velocities  behind  moving  vessels. 

Method  1:  Fuehrer.  Romisch.  and  Eneelke 

7.  Equations  are  presented  by  Fuehrer,  Romisch,  and  Engelke  (1981)  for 
estimating  the  maximum  bottom  velocity  for  moving  vessels.  Power  functions 
are  used  instead  of  the  exponential  forms  and  are  given  as 


Vb(x) 

V0 


(6) 


where 

Vb(x)  =  bottom  velocity  in  x-direction  at  coordinate  x 
V0  =  propeller  jet  velocity  at  x  =  0 
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A  —  function  of  propeller  height  above  bottom  and  if  rudder  is 
present  behind  propeller 

Dp  =  propeller  diameter 

a  -  0.6  if  spreading  is  limited  by  bottom  and  water  surface 
-  0.3  if  spreading  is  also  limited  by  adjacent  wall 

The  propeller  jet  velocity  can  be  computed  from 

Vc  r  1.6nDp(K,o>0-5 


where 

n  -  propeller  speed,  revolutions  per  second 
Kto  —  thrust  coefficient  at  zero  speed  of  advance 
Often  Kto  is  unknown  and  Fuehrer  and  Romisch  (1977)  proposed  the  relation 


VQ  =  0 . 9  5nDp 


(8) 


to  obtain  VQ  to  within  +20  percent.  Blaauw  and  van  de  Kaa  (1978)  proposed 
the  following  relation,  which  is  not  dimensionless 


(9) 


where 

V0  =  propeller  jet  velocity,  m/sec 

C'  =  1.17  for  ducted  propellers 
=  1.48  for  nonducted  propellers 

P  =  engine  power,  kW 

Dp  =  propeller  diameter,  m 

Fuehrer,  Romisch,  and  Engelke  (1981)  introduce  a  modified  advance  coefficient 
J  defined  as 


J  = 


V 


nD„ 


(10) 
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Using  the  modified  advance  coefficient,  the  relation  for  maximum  bottom  veloc¬ 
ity  for  maneuvering  versus  moving  navigation  was  found  to  be 


^b.max  ^b,  max  ,  J*0  (  1  ”3) 


(ID 


where 

Vbmax  =  maximum  bottom  velocity  for  moving  navigation 

Vb  maXjJ=o  =  maximum  bottom  velocity  for  maneuvering  navigation 
(vessel  speed  =  0) 

The  maximum  bottom  velocity  for  maneuvering  navigation  is 


V, 


b.max, J=0 


EV0 


wr 


where 

E  =  0.25  for  inland  vessel,  tunnel  stern,  twin  rudder  gear 
hp  =  distance  from  center  line  of  propeller  to  bottom 
Combining  Equations  11  and  12  leads  to 


V  =  F 

vb,majc 


D„ 


-1.0 


v„(i  -  J) 


(12) 


(13) 


Note  that  this  technique  addresses  only  the  maximum  bottom  •  elocity,  not  the 
velocity  field,  which  is  of  interest  in  this  study.  This  technique  does  dem¬ 
onstrate  the  reduction  in  velocity  that  occurs  with  increasing  vessel  speed. 
Method  2:  Verhev 

8.  Verhey  (1983)  developed  a  method  for  determining  the  propeller 
velocities  behind  a  moving  vessel.  The  induced  jet  velocity  for  a  moving 
vessel  is 


Vo  =  K  +  Kt 


8nzDi 


0.5 


(1^) 
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where 


Kt  =  thrust  coefficient  (depends  on  ship  speed) 
n  -  3.1416 

Va  =  entrance  velocity  defined  as 


va  =  (V  ♦  Vr)(l  -  W) 


(15) 


where 


Vr  =  return  velocity  as  determined  from  Schijf  (1949), 
Bouwmeester  et  al .  (1977),  or  other  methods 

W  =  wake  fraction  (0.3-0. 5  for  push  tows) 

The  contraction  diameter  D0  is 


D 


2 

o 


D2 

K.4n22 

7T 

v0<v0  -  Va) 


(16) 


9.  The  next  step  is  determining  the  wake  velocity,  which  is  the  veloc¬ 
ity  of  the  water  set  in  motion  behind  the  vessel  acting  in  the  same  direction 
as  the  vessel  is  traveling.  The  wake  velocity  is  assumed  constant  below  a 
horizontal  plane  at  Z/D0  =0.5  ,  where  Z  is  the  vertical  coordinate  mea¬ 
sured  from  the  center  line  of  the  propeller.  Determination  of  the  wake  flow 
is  based  on  the  diffusion  theory  for  jets  given  by  Albertson  et  al .  (1950). 
The  wake  flow  equations  are  as  follows: 

a.  The  zone  of  flow  establishment  (x/D0  <  6.2): 


Vw  «  V  -  (V  -  Va)  exp 


f_L  +  0.081 21  -  0.5]2 

Dn  D„ 


(17) 


where  Vw  is  the  velocity  of  wake  below  Z/Dc  =  0.5 
b.  The  zone  of  established  flow  (x/D0  >  6.2): 
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10. 

desired 


'  r  ' 

2  ' 

v.  ■  v  -  O'  -  exp 

-76.2  - 

K 

X 

1 

% 

(18) 


with 


(V  -  VJ 

'  a/ma* 

"Tv-3  va)“ 


=  6.2 


-1.0 


(19) 


The  next  step  is  to  determine  the  propeller  jet  velocity  at  the 
x  and  r  : 

a.  The  zone  of  flow  establishment: 

Let  the  velocity  increment  U2  be  defined  as 


U2  =  VQ  -  Va  (20) 

Jet  velocity  is 


where 

r-b.  <)"' 

y  =  horizontal  distance  from  the  center  line  of  the 
propeller 

C  -  0.18 

b.  The  zone  of  established  flow: 


V  (x). 


(22) 
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where  V(x)max  is  the  maximum  velocity  at  specified  x  ,  r  =  0 
defined  as 


V(x).4X 

u2 


(23) 


where 

b'  -  1  +  J. 

Ja  -  advance  coefficient  defined  as 


J 


(24) 


The  resulting  bottom  velocity  in  the  propeller  wash  region  Ve  is 


V,  =  vx 


+  v„  -  V 


(25) 


According  to  Verhey  (1983) 

It  must  be  stated  that  the  method  presented  for  calcu¬ 
lation  of  the  velocities  in  the  propeller  jet  behind  a 
sailing  vessel  is  rather  rough,  compared  to  the  method 
used  to  calculate  the  velocities  behind  maneuvering  ves¬ 
sels.  The  uncertainties  in  the  wake  fraction  W  ,  the 
impossibility  of  measuring  the  velocities  induced  by  the 
propeller  alone  and  the  schematizing  of  the  flow  field 
are  some  of  the  problems  to  be  solved.  Hence,  the  method 
presented  can  only  give  an  indication  of  the  velocities 
to  be  expected.  It  will  be  obvious  that  the  velocities 
behind  a  maneuvering  ship  will  be  greater  in  all  cases. 

Method  3:  Oebius 

11.  Oebius  (1984)  presents  a  method  for  determining  the  velocity  field 
behind  a  moving  vessel.  The  basic  equations  are  as  follows: 
a.  The  diameter  Dc  defined  as: 


0o  =  2(0 . 67Rf  «■  RJ 


(26) 
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where 


R{  -  blade  radius  (from  outside  of  hub  to  blade  tip) 
—  hub  radius 

b.  The  zone  of  flow  establishment  (x  <  xc)  : 


Vx  =  Vc  exp 


-0.5 

'  r  -  DX  ' 

2 

-  DX 

»  j 

(27) 


with  the  limit  of  the  flow  establishment  zone  xD  given  by 


=  2Do  exP 


1.265 


(28) 


where 

D„  =  modified  propeller  diameter 

Vi  =  velocity  at  infinity  assumed  equal  to  vessel  speed 

DX  =  distance  from  propeller  center  line  to  maximum 
velocity  given  by 


DX  =  0.32  Dq 


ur 


(29) 


and  the  following  equation  which  is  not  dimensionally  correct 


d;  =  1 . 15D„V. 


0.055 


(30) 


c.  The  zone  of  established  flow  (x  >  xc)  : 


Vx  =  V(x)max  eXP 


(31) 
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with 


VCx)^ 


with 


/3  =  0 . 6  exp 


and 


X  = 


K 

T 


-  DX  +  0.0875(x  -  xD) 


(32) 


(33) 


(34) 


and  xQ  according  to  Equation  28. 


Method  4:  Balanin  and  Bvkov 

12.  Balanin  and  Bykov  (1965)  presented  a  method  based  on  a  free  turbu¬ 
lent  jet.  However,  one  of  the  terms  was  not  defined  and  this  method  was  not 
pursued  in  this  study. 

Method  5:  Hochstein 

13.  Hochstein  and  Adams  (1986)  present  an  equation  for  propeller  jet 
velocities  derived  from  equations  presented  by  Blaauw  and  van  de  Kaa  (1978). 
The  Blaauw  and  van  de  Kaa  (1978)  equations  were  developed  for  vessel  speed 
close  to  zero  but  Hochstein  and  others  have  used  the  equation  in  a  wide  vari¬ 
ety  of  studies  for  vessels  fully  underway.  The  basic  equation  is 


CNP 

dpP 

1/3 

r 

(35) 


where 

V(r)  =  maximum  propeller  jet  velocity  at  radial  distance  r  from 
propeller  axis,  fps 
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C  «  coefficient  =  1.285 
Np  =  number  of  propellers 

P  «  total  engine  power,  hp 

The  difficulty  in  using  this  equation  for  tows  underway  is  that  the  coeffi¬ 
cient  C  must  be  determined  for  moving  vessels.  USAED,  Huntington  (1980b), 
references  measurements  on  the  Ohio  River  by  the  Huntington  District  in  the 
determination  of  C  -  1.285  .  However,  most  of  the  measurements  were  far  from 
the  towboat  and  probably  were  not  the  result  of  the  propeller  jet.  Hochstein 
and  Adams  (1986)  reference  measurements  on  the  Kanawha  River  that  verify 
C  =  1.285  .  Tow  speed  is  implied  in  the  coefficient  C  .  The  advantage  of 
this  method  lies  in  its  simplicity  and  ease  of  application. 

Displacement  and  Return  Velocity  Studies 

14.  In  addition  to  velocities  induced  by  the  propeller  jet,  the  dis¬ 
placement  of  water  by  the  moving  vessel  can  also  create  significant  vessel- 
induced  velocities  whose  magnitude  is  primarily  dependent  on  vessel  speed, 
average  channel  depth,  and  the  blockage  ratio.  The  primary  displacement- 
induced  velocity  acts  opposite  to  the  direction  of  travel  and  is  referred  to 
as  "return  velocity."  The  displacement-induced  return  velocity  is  also  accom¬ 
panied  by  a  lowering  of  the  water  level  between  vessel  and  bank,  which  is 
referred  to  as  "drawdown."  This  drawdown  is  greatest  near  the  vessel  and  is 
responsible  for  vessel  squat,  which  causes  a  vessel  to  have  a  reduced  under¬ 
keel  clearance  when  it  is  underway.  Numerous  techniques  are  available  for 
determining  the  average  return  velocity  and  average  water-level  drawdown. 

Three  of  these  will  be  presented  in  the  following  paragraphs.  As  given  in 
USAED,  Huntington  (1980a),  the  approaches  are  one-dimensional  and  it  is  neces¬ 
sary  that  certain  assumptions  be  made: 

a.  Constant  ship  speed,  in  a  channel  of  uniform  trapezoidal  or 
rectangular  cross  section. 

b.  Straight  channel  of  infinite  length. 

c.  Uniform  cross  section  of  ship,  disregarding  shape. 

d.  Uniform  return-current  velocity  around  ship  in  channel  cross 
section . 

e.  Uniform  water-level  depression  alongside  ship  in  channel  cross 
section . 

f.  Squat  over  ship's  length  equal  to  water-level  depression. 
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g.  Friction  losses  disregarded. 

The  Permanent  International  Association  of  Navigation  Congresses  (PIANC) 
(1987)  states  that  most  methods  are  limited  to  waterway  width/beam  width 
ratios  of  2-12  and  recommended  Schijf  for  loaded  pushtows  and  Bouwmeester  for 
other  ship  types. 

Method  1:  Schiif 

15.  Schijf  (1949)  used  a  conservation  of  energy  approach  and  developed 
the  equation  for  average  water-level  drawdown  as 


1  - 


1 

N 


z 

’H 


■ 

Y1 

Z 

1*0 

T5 

1  +  Z 

V4 

WJ 

(36) 


where 

z  =  average  water-level  drawdown 
h  =  average  channel  depth  -  area/top  width 
g  =  acceleration  due  to  gravity 
and  the  equation  for  average  return  velocity  Vt  as 


1  - 


1 

N 


V2 

2gh 


(37) 


The  equations  have  been  graphically  solved  and  are  presented  in  Jansen  and 
Schijf  (1953).  PIANC  (1987)  presents  an  alpha  factor  used  to  correct  the 
Schijf  method  to  improve  the  comparison  between  observed  and  computed  values 
of  drawdown  and  return  velocity.  At  high  enough  vessel  speeds,  return 
velocities  become  large  enough  to  reach  critical  conditions  and  a  self- 
propelled  vessel  reaches  its  so-called  limiting  velocity,  which  cannot  be 
exceeded.  The  Jansen  and  Schijf  (1953)  equation  for  limiting  velocity  is 


1  1  vc2t 

"  TJ  +  7  "gK 


3 

7 


11/3 


(38) 
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where  Vcr  is  the  limiting  velocity  for  self-propelled  ships.  This  relation 
has  been  verified  in  canals  having  relatively  small  values  of  N  (<10)  . 
Method  2:  Bouwmeester 

16.  Blaauw  and  van  der  Knaap  (1983)  provide  a  comprehensive  evaluation 
of  the  various  methods  and  conclude  that  the  Bouwmeester  et  al .  (1977)  rela¬ 
tion  provides  the  best  estimates  for  average  return  velocity  and  average 
water-level  drawdown.  The  Bouwmeester  relation  is  based  on  a  conservation  of 
momentum  and  water-level  drawdown  as  shown: 

0 . 5 

(39) 
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B  =  0.3  e  u-0(v/vcr,)  if  V/Vcr  <  0.65 
e  -  2.7183 

B  -  1  if  0.65  <  V/Vct  <  1 

For  Ohio  River  studies,  B  =  1  was  used  by  Louis  Berger  and  Associates 
(USAED,  Huntington,  1980b),  and  the  basic  equation  becomes 


vr  =  v  (yr  - 1) 

The  limiting  speed  VCI  is  determined  by  Hochstein  as 


(92) 


Vcr  =  K 


0.5 


(93) 


where 

K  =  constrainment  factor  given  in  Hochstein  and  Cohen  (1980) 

A0  =  cross-section  area  of  the  waterway 
18.  The  same  method  for  determining  the  limiting  velocity  VCI  is  given 
in  Fuehrer  and  Romisch  (1977) .  A  graph  is  provided  giving  the  constrainment 
factor  K  as  a  function  of  (h/d) (L/b)  and  (L/B0)  ,  where  L  is  the  vessel 
length.  The  Jansen  and  Schi j f  (1953)  and  Fuehrer  and  Romisch  (1977)  relations 
for  limiting  velocity  are  compared  in  the  following  tabulation  for  conditions 
similar  to  those  on  the  Ohio  River: 


L 

b 

h 

d 

B 

O 

V 

cr 

ft* 

ft 

ft 

ft 

ft 

N 

fos 

1,000 

105 

20 

9 

1,800 

38.1 

21** 

20  f 


*  A  table  of  factors  for  converting  non-SI  units  of 
measurement  to  SI  (metric)  units  is  found  on  page  3. 
**  Jansen  and  Schijf  (1953) 
f  Fuehrer  and  Romisch  (1977) 


For  a  loaded  tow  sailing  on  the  channel  center  line,  both  methods  give  similar 
results.  The  Fuehrer  and  Romisch  (1977)  study  also  describes  tests  to  deter¬ 
mine  Vcr  for  vessels  sailing  off  the  center  line  in  a  uniform  channel: 
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Vcr (of f  axis)  =  VCI(1  -  0.15a'  ) 


(44) 


where 

a'  -  eccentricity  factor  for  vessel  sailing  off  channel  center  line, 
zl/(0.5Bo) 

zl  —  distance  from  canal  axis  to  vessel 
Fuehrer  and  Romisch  also  developed  equations  for  the  displacement  velocity 
beneath  the  vessel 


a 


o 


..11  +  5.25 


1 

2.2 

TJ 
.  > 

(45) 


where 

actual  discharge  under  vessel  bottom 

0  ~  - vFE - 

(V  -  VM)  [h  -  (d  +  z)  ]b 

"  vHB - 


where  Vbd  is  the  displacement  velocity  beneath  the  vessel.  Note  that  in 
Equation  45  the  influence  of  channel  size  is  small  except  for  very  low 
blockage  ratios. 

19.  The  variation  of  the  return  velocity  from  vessel  to  bank  is  also  of 
importance  in  this  study.  Fuehrer  and  Romisch  (1977)  present  an  equation  for 
determining  the  maximum  return  velocity 


a 


(46) 


with 


a 


max 


1,0.114 


B0 

TT 


0.715 


(47) 
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where  Vr  max  is  the  maximum  return  velocity  at  the  vessel,  a  is  an  empir¬ 
ical  shape  factor  used  to  increase  the  average  return  velocity  obtained  from 
either  the  Bouwmeester,  Schijf,  Hochstein,  or  other  return  velocity  relations. 
For  a  <  1.5  ,  the  shape  of  the  return  velocity  distribution  should  be  repre¬ 
sented  by  a  linear  function  (USAED,  Huntington  1980b).  Louis  Berger  and  Asso¬ 
ciates  (USAED,  Huntington,  1980b)  uses  the  following  equation  to  determine  the 
distribution  of  return  velocity  from  vessel  to  bank  for  a  >  1.5 


Vr(y)  =  Vr iMX  exp  (-y/K2)  (48) 

where 

Vr(y)  =  return  velocity  as  a  function  of  distance  from  vessel 
K2  =  Bside  /  (a  {  1  -  exp  [  -F  (a)  a] }) 

Bside  =  distance  from  vessel  to  bank  (must  be  >B0/6  according  to 
USAED,  Huntington  (1980b)) 

F(a)  =  0.42  +  0.52  In  a 

20.  McNown  (1976)  states  that  as  the  ratio  N  of  waterway  cross-section 
area/vessel  cross-section  area  becomes  large  (not  quantified) ,  the  distribu¬ 
tion  of  return  velocity  becomes  nonuniform.  Delft  Hydraulics  Laboratory 
(1979)  presents  the  distribution  of  return  velocities  for  a  vessel  moving  both 
on  the  channel  center  line  and  close  to  one  bank.  Results  showed  a  maximum 
velocity  beneath  and  close  beside  the  vessel.  The  maximum  return  velocity  was 
located  near  the  bow  for  the  vessel  on  the  center  line  and  near  the  rear 
one-third  point  for  the  vessel  near  the  bank.  These  results  agree  with  the 
results  of  McNown  (1976)  regarding  uniformity  of  the  velocity  distribution  for 
vessels  near  the  bank  and  nonuniformity  of  velocity  distribution  for  vessels 
far  away  from  the  bank. 

21.  Most  of  the  equations  for  return  velocity  assume  that  the  vessel  is 
sailing  on  the  channel  center  line.  Of  importance  to  this  study  are  also  con¬ 
ditions  where  vessels  sail  off  the  channel  center  line.  Marchal  and  Spronck 
(1977)  determine  the  return  velocities  for  this  condition  by  treating  each 
side  of  the  vessel  as  being  independent  of  the  other  side.  To  determine 
return  velocities,  the  waterway  area  used  in  the  predictive  equations  is  equal 
to  two  times  the  area  between  the  vessel  center  line  and  the  bank  on  the  side 
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for  which  the  return  velocities  are  being  determined.  For  the  side  where  the 
vessel  is  close  to  the  bank,  the  waterway  area  used  in  the  return  velocity 
equations  will  be  less  than  the  actual  waterway  area.  For  the  side  with  the 
greater  distance  from  vessel  to  bank,  the  waterway  area  used  in  the  equations 
will  be  larger  than  the  actual  waterway  area.  PIANC  (1987)  presents  an  equa¬ 
tion  for  determining  the  effective  area  for  vessels  sailing  off  the  channel 
center  line  in  a  prismatic  channel.  Blaauw  et  al.  (1984)  present  a  plot  show¬ 
ing  the  variation  of  maximum  return  velocity  with  eccentricity. 

22.  The  variation  of  the  return  velocities  with  distance  from  the  vessel 
may  be  accompanied  by  a  variation  in  the  vertical  velocity  distribution. 
However,  Fuehrer  and  Romisch  (1977)  observed  a  nearly  uniform  velocity 
distribution  from  bottom  to  water  surface  in  the  return  velocity.  They  also 
observed  a  transition  region  next  to  the  vessel  that  had  transverse  veloc¬ 
ities.  Blaauw  et  al .  (1984)  used  the  Schlichting  (1968)  formula  for  flow  over 
rough  plates  to  define  the  relationship  between  return  velocity  and  shear 
stress . 

Tow-Induced  Velocities  Measured  in  Field  Studies 

23.  USAED ,  Huntington  (1980a),  conducted  measurements  of  about  200  tows 
on  the  Ohio  River  as  part  of  a  study  concerning  the  replacement  lock  at  Galli- 
polis.  Two-dimensional  electromagnetic  meters  were  used  in  shallow-water 
areas  near  the  bank,  and  a  three-dimensional  electromagnetic  meter  was  used  on 
the  river  bottom  near  the  path  of  the  vessels.  (One  tow  came  within  16  ft  of 
the  meter.  All  other  tows  were  52  ft  or  greater  from  the  three-dimensional 
meter.)  The  Huntington  study  recommends  either  Hochstein's  or  Bouwmeester ' s 
equation  for  return  velocity.  The  supplement  to  the  USAED,  Huntington 
(1980a),  report  by  Louis  Berger  and  Associates  (USAED,  Huntington,  1980b) 
states  that  similar  results  are  obtained  from  Hochstein  and  Schijf  equations 
for  return  velocity.  The  Hochstein  relation  for  return  velocity  was  multi¬ 
plied  by  a  factor  of  1.1  to  obtain  better  agreement  with  measured  values.  The 
USAED,  Huntington  (1980a),  report  recommends  the  Blaauw  and  van  de  Kaa  (1978) 
equation  for  propeller  jet  velocity.  This  report  states  that  the  area  of 
impact  from  propeller  jet  velocities  is  60-65  ft  wide.  Open-wheel  propulsion 
systems  were  found  to  increase  turbidity  more  than  Kort  nozzle  systems. 

24.  Environmental  Science  and  Engineering  (ESE)  (1981)  conducted  lield 


measurement  of  velocities  for  30  tows  on  the  upper  Mississippi  River  and 
29  tows  on  the  Illinois  River.  Two-dimensional  electromagnetic  meters  were 
installed  1  ft  above  the  bottom  at  two  positions  adjacent  to  the  navigation 
channel,  one  near  the  shore  and  the  other  near  the  sailing  line.  ESE  reported 
that  the  effects  of  a  tow  were  generally  measurable  more  than  4  min  before  the 
tow  reached  the  instruments.  ESE  reported  that  the  Hochstein  relation  for 
return  velocities  underestimated  the  measured  bottom  velocities  by  a  factor  of 
2  for  nearshore  velocities  and  by  an  average  of  30  percent  for  offshore  veloc¬ 
ities.  A  correction  factor  was  developed  to  improve  the  measured  versus 
computed  results.  ESE  pointed  out  that  all  comparisons  were  based  on  bottom 
velocities  and  that  depth-averaged  velocities  would  also  be  underestimated. 

ESE  reported  on  a  comparison  of  velocities  at  two  depths  (total  depth  14.5  ft) 
that  showed  that  the  surface  return  velocity  (6  ft  belov  surface)  due  to  tow 
passage  was  1.5  times  the  bottom  velocity  (1  ft  above  bottom)  due  to  tow 
passage . 

25.  Bhowmik  (1981)  conducted  field  measurements  of  velocity  for  19  tows 
on  the  uoper  Mississippi  and  22  tows  on  th  *inois  Rivers.  One-dimensional 
Price  current  meters  were  used  to  obtain  velocities  at  depths  of  0.95,  0.8, 
0.6,  and  0.2  of  the  depth  measured  from  the  surface. 

26.  Hochstein  and  Auams  (10°4)  reported  on  velocity  measurements  taken 
on  the  Kanawha  River  to  check  the  Hochstein  relations  for  return  and  propeller 
jet  velocities.  Results  showed  good  agreement  between  the  observed  data  and 
the  Hochstein  relations. 

Recommendations  for  Study 


27.  Based  on  analysis  of  the  existing  literature,  the  following  list 
recommends  areas  in  which  additional  study  is  needed: 

a.  Relations  for  propeller  jet  velocity  for  moving  tows  should 
consider  speed  in  the  analysis.  Verhey  (1983)  is  the  most  com¬ 
prehensive  relation  but  also  requires  the  most  input.  The  width 
of  the  propeller  jet  attack  on  the  bottom  for  a  moving  tow  needs 
to  be  defined. 

b.  Existing  return  velocity  relations  do  not  provide  satisfactory 
comparisons  with  field  data  gathered  in  large  rivers.  Return 
velocity  distribution  equations  need  to  be  developed  for  tows 
moving  in  asymmetric  cross  sections  at  various  positions  across 
the  section.  The  influence  of  ambient  currents  and  tow  direc¬ 
tion  needs  to  be  evaluated. 
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c.  The  use  of  velocity  in  studies  of  navigation  effects  needs  to  be 
standardized.  Consideration  should  be  given  to  how  the  veloc¬ 
ities  are  going  to  be  used  in  the  navigation  effects  studies. 

For  example,  in  studies  of  sediment  movement  induced  by  naviga¬ 
tion,  most  existing  sediment  transport  relations  are  based  on 
either  shear  stress  or  depth-averaged  velocity.  Depth-averaged 
velocity  is  valid  for  return  velocities  in  the  area  between  the 
vessel  and  the  bank.  However,  depth-averaged  velocity  is  mean¬ 
ingless  in  areas  such  as  the  propeller  jet  because  bottom  and 
surface  velocities  may  be  in  different  directions.  In  this 
case,  techniques  are  needed  to  transfer  bottom  velocity  to 
either  shear  stress  or  an  equivalent  depth-averaged  velocity. 


24 


PART  III:  FLOW  VISUALIZATION  STUDIES  WITH  1: 20-SCALE  MODEL 

28.  Flow  visualization  studies  were  conducted  to  gain  a  better  under¬ 
standing  of  the  flow  patterns  near  a  moving  vessel.  The  flume  used  in  these 
studies  represented  slack-water  conditions  having  zero  ambient  flow.  Results 
from  this  study  were  primarily  in  the  form  of  video  tapes  taken  both  under¬ 
water  and  from  above  the  moving  tow.  Velocity  vector  plots  for  two  of  the 
test  conditions  and  movement  tests  with  lightweight  plastic  beads  were  also 
determined. 

29.  The  physical  model  used  in  this  study  was  constructed  to  a  scale 
ratio  of  1:20.  The  model  layout  is  shown  in  Plate  1.  A  winch  was  located  at 
one  end  of  the  model  for  use  in  the  bead  tests  only.  An  array  of  nails  was 
installed  on  the  floor  of  the  model  and  yarn  strings  approximately  0.25  ft 
(model)  in  length  were  attached  to  the  nails  2  ft  above  the  bottom  (all  dimen¬ 
sions  are  in  prototoype  unless  stated  otherwise) .  The  yarn  strings  had  a 
fairly  significant  resistance  to  movement,  and  other  flow  visualization  tech¬ 
niques  were  tried.  The  best  techniques  found  were  dye  and  lightweight  plastic 
beads.  The  dye  was  injected  through  1/4-in. -diam  (model)  copper  tubing  placed 
on  the  channel  bottom  with  1/64-in. -diam  (model)  holes  drilled  in  the  top  of 
the  copper  tubing.  The  plastic  beads  were  used  in  qualitative  movement  tests 
that  will  be  discussed  in  a  later  paragraph.  The  bottom  of  the  model  was 
painted  white,  and  underwater  lighting  was  used  to  improve  light  conditions 
beneath  the  tow.  A  cable  was  placed  the  full  length  of  the  model  flume  along 
the  sailing  line  of  the  tow,  and  guides  on  the  bow  of  the  lead  barge  and  the 
stern  of  the  towboat  were  used  to  ensure  a  consistent  sailing  line.  Depth  of 
water  ranged  from  15  to  30  ft. 

30.  The  1: 20-scale  towboat  represented  a  5,600-hp  towboat  with  twin 
9-ft-diam  propellers,  Kort  nozzles,  and  main  and  flanking  rudders  in  line  with 
each  propeller  shaft.  The  towboat  dimensions  were  45.6  ft  wide  by  209  ft 
long,  and  the  towboat  draft  was  9  ft  for  all  tests.  All  tests  were  conducted 
with  a  0-deg  rudder  setting,  and  both  propellers  were  turning  the  same  speed. 
Looking  at  the  stern  of  the  towboat,  the  starboard  propeller  was  turning  coun¬ 
terclockwise  and  the  port  propeller  was  turning  clockwise.  Propeller  speeds 
could  be  varied  from  130  to  190  rpm.  The  35-ft-wide  by  195-ft-long  barges 
drafted  9  ft  when  leaded  and  1 . 5-2 . 0  ft  when  unloaded.  The  individual  barges 
were  lashed  together  to  form  a  tow  up  to  three  wide  by  three  long.  The  bows 
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of  the  lead  barges  were  raked  on  a  radius  of  25  ft.  The  25-ft  radius  was 
extended  from  the  bottom  of  the  barge  for  a  vertical  distance  of  10  ft.  The 
stern  of  the  rear  barges  had  boxed  ends.  The  model  is  shown  in  Figure  1. 

31.  To  ensure  the  similarity  of  flow  patterns  between  model  and  proto¬ 
type,  the  model  was  operated  with  an  equal  Froude  number  in  model  and  proto¬ 
type.  The  following  relations  were  used  to  transfer  quantities  from  model  to 


prototype : 

Characteristic 

Dimension* 

Scale  Relations 

Model : Prototvoe 

Length 

1:20.0 

Weight  or  volume 

4 

1:8,000.0 

Time 

U/2 

1:4.4721 

Velocity 

4/2 

1:4.4721 

rpm 

i/4/2 

1:0.224 

*  Dimensions  are  in 

terms 

of  length  ratio. 

Reduced  scale  navigation  models  have  proportionally  greater  frictional  forces 

than  in  the  prototype.  This  results  in  slower  tow  speeds  in  the  model  than  in 

the  prototype.  For  these  flow  visualization  studies,  the  speed  difference  was 

not  considered  to  affect  results.  Several  tests  were  conducted  with  a  towing 

mechanism  providing  all  the  propulsion. 


32.  The  following  observations  were  made  during  the  flow  visualization 
tests : 

a.  In  front  of  the  tow,  water  was  set  in  motion  in  the  same  direc¬ 
tion  as  the  tow.  Just  in  front  of  the  bow  of  the  tow,  veloc¬ 
ities  at  the  water  surface  were  approximately  equal  to  the  tow 
speed.  At  the  channel  bottom,  these  velocities  were  reduced  and 
depended  on  the  depth  of  flow.  Velocities  2  ft  above  the  bottom 
at  the  bow  of  the  towboat  will  be  referred  to  as  Vbb  . 

b.  Just  downstream  of  the  bow  (beneath  the  tow) ,  there  was  a  rapid 
reversal  in  flow  that  resulted  in  flow  opposite  to  the  direction 
of  the  tow.  This  flow  is  related  to  the  displacement  effects  of 
the  tow,  and  velocities  will  be  referred  to  as  Vbd  .  The  maxi¬ 
mum  Vbd  occurred  approximately  40  ft  behind  the  bow  of  the 
lead  barge . 

c.  Vbd  decreased  toward  the  stern  of  the  rear  barge  for  the  15-ft 
depth  and  stayed  about  the  same  for  the  30-ft  depth.  Vbd  may 
be  influenced  by  the  suction  effects  of  the  propellers.  The 
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Figure 


magnitude  of  the  suction  effect  will  depend  on  several  factors 
including  the  length  of  the  towboat.  The  model  towboat  simu¬ 
lates  a  relatively  long  towboat  (209  ft)  and  suction  effects 
will  be  less  than  for  the  typical  5,600-hp  towboat,  which  is 
about  150  ft  in  length. 

d.  From  the  bow  to  the  stern  of  the  towboat,  there  existed  a  highly 
complex  flow  field,  based  on  observation  of  the  dye  movement  be¬ 
neath  the  towboat.  Velocities  in  this  region  are  influenced  by 
the  wake  from  the  upstream  barges,  suction  effects  of  the  pro¬ 
peller,  tow  speed,  depth,  and  other  factors.  At  the  15-ft 
depth,  flow  patterns  did  not  indicate  significant  velocities  in 
a  horizontal  plane  from  the  towboat  bow  to  the  propellers.  How¬ 
ever,  the  dye  pointed  up  as  the  propellers  passed  over  the  dye 
location.  Any  velocity  measurements  taken  in  this  region  must 
be  three-dimensional  to  be  valid. 

e.  Behind  the  propellers,  the  dye  pointed  in  the  direction  of  the 
tow  prior  to  the  arrival  of  the  propeller  jet.  This  flow,  which 
is  in  the  same  direction  as  the  tow,  results  from  infilling  be¬ 
hind  the  towboat  and  is  referred  to  as  a  wake  flow. 

f.  In  the  propeller  wash  region,  the  velocity  field  is  also  complex 
and  is  generally  opposite  the  direction  of  tow  travel.  Veloci¬ 
ties  are  influenced  by  the  wake  of  the  barges  and  towboat,  pro¬ 
peller  jet,  tow  speed,  depth,  and  other  factors.  Bottom  veloci¬ 
ties  in  this  zone  will  be  referred  to  as  V„  .  The  width  of  the 
propeller  jet  attacking  the  bottom  was  observed  to  be  about 

50  ft,  which  compares  well  with  the  60-65  ft  reported  in  the 
Gallipolis  study  (USAED,  Huntington,  1980a) .  Outside  this 
region,  the  wake  flow  behind  the  barges  of  the  tow  creates  a 
velocity  field  in  the  same  direction  as  the  tow. 

g.  Several  tests  were  conducted  with  the  propellers  not  spinning 
and  the  towing  mechanism  providing  all  the  propulsion.  These 
tests  demonstrate  that  the  tow  sets  a  large  volume  of  water  into 
motion  behind  the  tow  and  in  the  same  direction  as  the  tow. 

This  water  motion  is  opposite  to  the  direction  of  the  propeller 
jet  and  becomes  significant  for  three-wide  loaded  tows  in 
shallow  water.  These  effects  are  accounted  for  by  the  term  Vw 
in  the  Verhey  (1983)  method. 

33.  Velocity  vector  plots  were  prepared  for  tow  configurations  of 
three  wide  by  three  long  loaded  (Plate  2)  and  three  wide  by  three  long  un¬ 
loaded  (Plate  3).  The  observations  described  in  paragraphs  32a-g  are  shown  in 
Plate  2.  These  represent  approximate  bottom  velocity  directions  and  magni¬ 
tudes  adjacent  to  the  moving  tow  based  on  observation  of  the  dye  patterns. 
Propeller  speed  for  both  vector  plots  was  185  rpm. 

34.  Lightweight  plastic  beads  having  a  specific  gravity  of  1.03  were 
used  to  conduct  qualitative  estimates  of  the  movement  of  sediment.  The  plas¬ 
tic  beads  had  a  fall  velocity  in  the  1 : 20-scale  model  equivalent  to  about  a  3- 
to  4-mm  quartz  sand  particle  in  the  prototype.  One  hundred  beads  were  placed 
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in  a  grouping  about  2  in.  (model)  in  diameter  at  various  locations  with 
respect  to  the  channel  center  line.  The  location  of  the  beads  was  recorded 
after  each  passage  of  the  tow.  The  following  conditions  were  tested  with  all 
tows  being  loaded: 


Configuration 
Width  (W)  by 
Length  (L) 

Tow  Speed 
mph 

Propeller 

Speed,  rpm 

Initial  Bead 
Location,  ft  Off 
Channel  Center  Line 

Plate 

3W  x 

3L 

6.9 

185 

12.5 

4 

3W  x 

3L 

6.9 

0  (towing  test) 

12.5 

5 

3W  x 

3L 

6.9 

185 

37.5 

6 

3W  x 

3L 

6.9 

185 

90.0 

7 

3W  x 

1L 

8.3 

185 

12.5 

8 

Most  of  the  plots  are  relatively  straightforward  except  for  the  beads  placed 
at  37.5  ft  off  the  channel  (Plate  6).  The  plot  indicates  that  the  beads  did 
not  move  significantly.  However,  the  underwater  camera  showed  that  the  beads 
moved  approximately  25  ft  opposite  to  the  tow  direction  due  to  the  displace¬ 
ment  flows  beneath  the  barges.  The  beads  were  then  redeposited  in  about  the 
original  location  by  the  wake  flow  behind  the  tow.  The  beads  at  37.5  ft  were 
not  niuved  by  the  propeller  jet. 
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PART  IV:  FLOW  VISUALIZATION  TESTS  OF  SURFACE  CURRENT 
PATTERNS  WITH  OLMSTED  MODEL 


35.  Flow  visualization  tests  were  conducted  in  the  upper  pool  of  the 
existing  1: 120-scale  navigation  model  of  the  Olmsted  Locks  and  Dam  on  the  Ohio 
River.  These  tests  were  conducted  to  determine  the  effectiveness  of  conven¬ 
tional  flow  pattern  visualization  techniques  as  applied  to  visualization  of 
flows  generated  by  moving  vessels .  These  tests  used  still  photography  with  a 
4-sec  (model)  shutter  speed  to  document  the  movement  of  confetti  during  pas¬ 
sage  of  the  tows  at  various  positions  in  the  channel.  The  overhead  camera  was 
positioned  16  ft  above  the  water  surface  of  the  model.  The  3/4-  by  3/4-in. 
(model)  confetti  was  placed  in  the  slack-water  pool,  and  the  test  was  con¬ 
ducted  only  after  all  movement  of  the  confetti  ceased.  Tests  of  two  channel 
sizes  were  conducted:  a  standard  15-barge  tow  (105  ft  wide  by  9-ft  draft  by 
1,150  ft  long)  operating  in  a  3 , 600-f t-wide  channel  and  a  standard  15-barge 
tow  operating  in  a  2 , 200— ft— wide  channel.  The  scale  ratios  were  1:120  and 
1:70  for  the  3 , 600-f t-wide  and  2 , 200-f t-wide  channels,  respectively.  Tests 
were  conducted  under  slack-water  conditions. 

36.  Model  quantities  were  transferred  to  the  prototype  by  means  of  the 
following  relations: 


Scale  Relations 
Model : Prototype 


Characteristic 

Dimension* 

1: 120-Scale  Model 

1:70— Scale  Model 

Length 

L  =  L  = 

R 

PrototVDe  Leneth 
Model  Length 

1:120 

1:70 

Velocity 

V  =  4'2 

1:10.95 

1:8.37 

Time 

T  =  4/2 

1:10.95 

1:8.37 

*  Dimensions  are  in  terms  of  length. 

37.  A  cross  section  and  tow  locations  for  the  1: 120-scale  channel  are 
shown  in  Plate  9.  Also  shown  are  the  limits  of  the  photographic  coverage, 
which  was  approximately  the  right  half  of  the  channel,  and  the  blockage  ratio 
N  .  The  following  tests  were  conducted: 
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Test 

Scale 

Ratio 

Tow 

Location 

Tow 

SDeed.  mDh 

Distance 
Moved,  ft 

Photo 

1 

1:70 

A 

7.6 

375 

1 

2 

1:70 

B 

7.5 

368 

2 

3 

1:70 

C 

7.0 

345 

3 

4 

1 : 120 

A 

10.0 

643 

4 

5 

1:120 

B 

8.7 

559 

5 

6 

1:120 

C 

7.6 

488 

6 

The  white  lines  perpendicular  to  the  tow  center  line  in  Photos  1-6  are  480  ft 
apart  in  the  1:120  scale  and  280  ft  apart  in  the  1:70  scale  and  can  be  used  to 
scale  distance  in  both  directions.  The  "distance  moved"  value  in  the  tabula¬ 
tion  refers  to  the  distance  the  tow  moved  during  the  4  sec  (model)  that  the 
shutter  of  the  camera  was  open.  For  example,  in  Photo  la,  the  0  shows  the 
position  of  the  bow  of  the  tow  when  the  shutter  was  opened.  The  C  in  Photo  lb 
shows  the  position  of  the  bow  when  the  shutter  was  closed.  Points  0  and  C 
were  375  ft  apart,  as  given  in  the  tabulation.  The  tow  traveled  a  significant 
portion  of  the  total  photograph,  which  is  important  to  remember  when 
evaluating  the  flow  patterns. 

38.  Analysis  of  the  photographs  shows  that  the  1: 70-scale  tow  produced 
surface  movement  from  tow  to  bank  even  for  position  A,  which  was  farthest  from 
the  bank.  Significant  movement  parallel  to  the  tow  resulting  from  return  ve¬ 
locities  was  observed  as  expected.  The  1 : 120-scale  tow  produced  surface  move¬ 
ment  only  in  a  width  of  about  300-400  ft  on  each  side  of  the  tow.  Surface 
movement  for  the  1: 120-scale  tow  was  predominantly  away  from  the  tow,  and  no 
significant  movement  parallel  to  the  tow  was  observed.  Quantitatively,  the 
photographs  for  the  1: 70-scale  tow  can  be  used  to  define  the  distribution  of 
return  velocities  at  the  surface.  Photos  1,  2,  and  3  were  used  to  define  the 
maximum  velocity  parallel  to  the  tow,  and  results  are  shown  in  Plates  10,  11, 
and  12,  respectively.  To  determine  velocity,  the  length  of  the  confetti 
streak  was  scaled  off  the  photograph.  This  length  was  divided  by  the  proto¬ 
type  time  that  the  camera  shutter  was  open  or  4(V ^0)  =  33.5  sec.  The  surface 
velocities  are  roughly  equal  to  the  depth-averaged  velocities  and  can  be  com¬ 
pared  to  the  computed  average  return  velocities  using  the  Bouwmeester  relation 

( SVeiTatTP  fnr  £»nt~irp  rrncc  Qprh'nn) 
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PART  V:  1: 20-SCALE  PHYSICAL  MODEL  INVESTIGATION  OF  BOTTOM 
VELOCITIES  NEAR  THE  PATH  OF  A  MOVING  TOW 

Description 

39.  The  1 : 20-scale  physical  model  tow  was  used  to  measure  velocities 
near  the  path  of  a  moving  tow.  A  sketch  showing  model  limits  is  shown  in 
Plate  13.  Also  shown  in  Plate  13  are  the  locations  of  the  velocity  meters 
used  in  the  study.  Three  channel  types  were  used  in  this  study  representing 
navigation  in  an  unconfined  channel,  navigation  near  one  bank,  and  navigation 
on  the  center  line  of  a  confined  channel.  A  channel  cross  section  at  the 
velocity  meter  location  shown  in  Plate  13  is  shown  in  Plate  14. 

40.  Details  of  the  1 : 20-scale  tow  are  given  in  paragraph  30.  The  thrust 

coefficient  for  the  propellers  at  zero  speed  of  advance  equals  0.51.  The  Kort 
nozzles  on  the  model  towboat  were  removed  to  conduct  several  open-  wheel  runs. 
Results  from  the  open -wheel  tests  were  qualitative  because  the  propellers  were 
designed  for  Kort  nozzles.  The  unconfined  channel  tests  were  conducted  with 
the  port  propeller  turning  counterclockwise  (as  viewed  from  the  rear  of  the 
towboat)  and  the  starboard  propeller  turning  clockwise.  At  this  point  in  the 
study  it  was  determined  that  the  propellers  of  most  towboats  rotate  in  the 
opposite  direction  to  that  used  in  the  unconfined  channel  tests.  The  propel¬ 
lers  were  switched  for  the  near-bank  and  confined  channel  tests. 

41.  The  scale  relations  are  the  same  as  those  given  in  paragraph  31. 

Because  excess  frictional  forces  are  present  in  scaled  navigation  models,  an 
added  force  is  required  in  the  model  to  obtain  equivalent  speeds  in  model  and 

prototype.  In  this  investigation,  a  towing  mechanism  was  used  to  provide  the 

added  force  and,  because  of  the  limited  model  length,  ensure  a  constant  vessel 
speed  while  the  model  tow  was  in  the  test  section. 

42.  Two-dimensional  electromagnetic  velocity  meters  were  initially 
placed  in  the  model,  but  the  measured  velocities  fluctuated  rapidly  when  the 
tow  was  in  the  vicinity  of  the  meters.  This  occurred  even  when  the  tow  was 
not  moving.  Further  investigation  of  the  velocity  meters  revealed  that  the 
frequency  response  was  not  fast  enough  to  monitor  the  rapid  changes  that  were 
taking  place  in  both  velocity  magnitude  and  direction.  A  Nixon  series 

400  propeller  meter  having  a  rotor  diameter  of  0.038  ft  (model)  was  tested  in 
the  model.  A  strip-chart  recorder  was  used  to  record  output  from  the  Nixon 
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meter.  Exhaustive  te'.,ts  to  determine  the  frequency  response  of  the  system 
were  not  conducted,  but  measured  velocities  were  compared  with  flow  patterns 
observed  beneath  the  vessel.  At  the  bow  of  the  loaded  tow  traveling  in  15  ft 
of  water,  a  rapid  flow  reversal  occurred  when  velocities  changed  from  positive 
to  negative  in  about  0.5  sec.  The  velocity  metering  system  captured  this  flow 
reversal,  showing  that  the  response  was  adequate  for  this  study.  The  propel¬ 
ler  meter  could  not  measure  direction  and  had  a  threshold  velocity  of  about 
0.5  fps  (0.1  fps  in  model). 

Test  Results 


43.  Tests  conducted  for  the  three  cross-section  types  are  summarized  in 
Tables  1-3.  The  limited  model  length  prevented  testing  tows  longer  than  two 
barges  long  plus  the  towboat.  In  many  tests  the  limited  model  length  required 
that  the  tow  be  stopped  and  the  test  ended  before  the  observed  velocites 
approached  zero.  The  observed  velocity  plots  shown  in  Appendix  A  are  ex¬ 
plained  by  a  master  l_g.'  .u,  and  the  relationship  between  velocity  meter  read¬ 
ing  and  prototype  \ e' .city  is  given  in  Table  Al .  Multiple  tests  were  con¬ 
ducted  for  each  combination  shown  in  Tables  1-3,  and  a  representative  test  was 
selected  for  each  combination.  The  observed  velocities  for  selected  tests  for 
the  three  cross-section  shapes  are  shown  in  Appendix  A.  Tests  1-40  are  for 
the  unconfined  channel.  Tests  114-231  are  for  the  near-bank  tests. 

Tests  300-330  are  for  the  confined  channel  tests. 

44.  A  series  of  tests  were  conducted  with  the  near-bank  channel  config¬ 
uration  to  determine  the  distance  from  the  tow  at  which  the  velocity  6  in. 
above  the  channel  bottom  was  equal  to  0.5  fps.  These  tests  were  conducted  on 
the  side  of  the  tow  opposite  the  near  bank.  Because  the  velocity  meter  could 
not  measure  velocities  this  low,  dye  was  injected  6  in.  above  the  channel  bot¬ 
tom  and  the  time  required  for  the  dye  to  traverse  a  fixed  distance  was  used  to 


determine 

these  velocities. 

Results  were 

as  follows: 

Test  No. 

DeDth.  ft 

Tow  Speed 
f  DS 

Distance  from 
Cable  or  Tow 
Center  Line,  ft 

Velocity  6  in. 
Above  Channel 
Bottom,  fps 

200 

15 

6.7 

300 

0.5 

201 

15 

11.2 

850 

0.5 

208 

30 

6.8 

130 

0.5 

209 

30 

11.4 

190 

0.5 
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45.  The  following  observations  were  made  about  the  model  velocities  near 
the  path  of  the  vessel: 

a.  For  equal  vessel  and  propeller  speeds,  the  unloaded  tow  produced 
higher  propeller  jet  bottom  velocities  than  did  the  loaded  tow. 
This  is  very  likely  the  result  of  the  wake  flow  (which  acts 
opposite  to  the  propeller  jet)  being  stronger  behind  the  loaded 
tow.  In  reality,  for  equal  propeller  speeds,  the  unloaded  tow 
would  be  going  faster  than  the  loaded  tow  and  propeller  jet 
velocity  at  the  bottom  for  the  unloaded  tow  would  be  reduced. 
This  reduction  would  be  the  result  of  two  factors.  (1)  the 
speed  of  the  propeller  jet  relative  to  the  bottom  would  be 
reduced;  and  (2)  the  higher  speed  would  produce  a  larger  wake 
flow,  also  counteracting  the  propeller  jet  velocity  at  the 
bottom. 

b.  For  loaded  tows,  the  maximum  center-line  displacement  velocities 
were  only  slightly  less  than  the  maximum  center-line  propeller 
jet  velocities  (both  2  ft  above  bottom).  However,  the  higher 
turbulence  intensity  of  the  propeller  jet  will  cause  the  pro¬ 
peller  jet  to  have  a  greater  transport  capacity. 

c.  At  equal  channel  depths,  bottom  center-line  displacement  veloc¬ 
ities  were  relatively  unaffected  by  the  change  in  channel  cross- 
section,  which  means  that  return  velocity  was  not  contributing 
to  the  total  velocity  for  the  blockage  ratios  used  in  this 
study . 

d.  Velocities  outside  the  vessel  (measured  27.5  ft  from  the  edge  of 
the  barge)  were  signif iciantly  influenced  by  blockage  ratio, 
which  means  that  return  velocity  was  contributing  to  the  total 
velocity. 

e.  The  towing  tests  (propeller  speed  of  0)  demonstrated  that  a 
large  wake  flow  develops  behind  a  vessel.  The  wake  flow  moved 
in  the  same  direction  as  the  vessel  and  persisted  for  relatively 
long  periods  of  time. 

f .  The  effects  of  tow  length  could  not  be  defined  because  of  the 
limited  range  of  this  parameter.  The  tests  of  open-wheel  versus 
Kort  nozzle  were  also  inconclusive  since  only  Kort  nozzle  type 
propellers  were  available. 

g.  The  direction  of  velocities  behind  the  towboat  measured  26  ft 
from  the  center  line  was  not  well  defined.  The  dye  injections 
demonstrated  that  this  is  a  borderline  area  for  the  influence  of 
the  propeller  jet,  which  means  that  the  observed  velocity  could 
have  been  a  wake  flow, 

h.  Bow  velocities  were  always  less  than  the  displacement  velocity, 
but  for  loaded  tows  traveling  in  shallow  water,  the  bow  veloc¬ 
ities  were  as  high  as  2.5  fps . 
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PART  VI:  PROTOTYPE  DATA  COLLECTED  BY  USAED,  LOUISVILLE 


46,  Prototype  tests  were  conducted  by  USAED,  Louisville,  in  1987  at  Ohio 
River  mile  581.  The  tests  were  conducted  with  leased  towboats  to  obtain  data 
regarding  velocities  near  the  river  bottom  at  a  limited  number  of  depths.  The 
details  of  the  prototype  tests  conducted  by  USAED,  Louisville,  are  shown  in 
Tables  4-6.  Three  different  tow  configurations  are  represented  by  the  proto¬ 
type  tests:  (a)  towboat  behind  one  unloaded  barge  ( Steve  Kuhr  (all  tests), 
John  Matthews  (all  tests),  Harold  Turner  (tests  01-04)),  (b)  towboat  behind 
one  loaded  barge  which  is  smaller  than  the  towboat  ( Harold  Turner 

(tests  05-08)),  and  (c)  towboat  operating  in  wake  region  of  three-wide  loaded 
barges  ( Harold  Turner  (tests  09-13)).  A  cross  section  showing  the  tow  loca¬ 
tions,  cross-sectional  areas,  and  widths  is  shown  in  Plate  15.  As  in  the 
physical  model,  measurement  of  velocities  beneath  a  moving  tow  in  the  proto¬ 
type  is  a  difficult  task.  Directions  and  magnitudes  change  rapidly  and  the 
two-dimensional  electromagnetic  velocity  meters  used  in  the  prototype  may  not 
have  been  fast  enough  in  their  frequency  response  to  capture  the  changes 
occurring  under  the  vessels.  The  meters  were  positioned  to  measure  velocities 
in  the  horizontal  plane  only. 

47.  A  summary  listing  of  the  observed  prototype  velocities  is  given  in 
Table  7.  This  listing  presents  the  velocity  magnitude  and  direction  for  five 
points  along  the  tow.  These  five  points  are  shown  in  Plate  2.  Each  entry 
represents  the  maximum  observed  for  that  location.  The  entries  with  a  ?? 
could  not  be  described  by  a  single  value.  The  following  presents  observations 
about  each  of  the  five  locations: 

a.  Bow  velocities .  Bow  velocities  are  those  ahead  of  a  tow.  For 
loaded  tows,  these  velocities  can  extend  a  considerable  distance 
in  front  of  the  tow.  For  the  four  tests  ( Harold  Turner 

tests  05-08)  having  a  or.e-wide  loaded  barge  (35  ft  wide)  in 
22  ft  of  water,  the  average  distance  ahead  of  the  tow  in  which 
the  bow  affected  the  measured  bottom  velocity  was  about  300  ft. 
For  the  five  tests  ( Harold  Turner  tests  09-13)  having  three- rfide 
loaded  barges  (105  ft  wide)  in  22  ft  of  water,  the  average  dis¬ 
tance  was  about  550  ft.  This  shows  that  for  the  tests  with  one 
unloaded  barge,  the  bow  velocity  caused  by  the  towboat  can  ex¬ 
tend  upstream  of  the  bow  of  the  unloaded  barge. 

b.  Displacement  velocity  near  bow  of  barge(s).  This  location  was 
selected  because  this  was  the  location  of  the  maximum  displace¬ 
ment  velocity  in  the  physical  model.  Similar  to  the  bow  veloc¬ 
ities,  the  displacement  velocities  for  the  tests  with  one  un¬ 
loaded  barge  were  affected  by  the  bow  effect  from  the  towboat. 
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Most  of  the  upbound  tows  caused  an  increase  in  the  bottom 
velocity  above  the  ambient  velocity  in  a  direction  opposite  the 
tow.  This  was  consistent  with  the  observation  of  displacement 
velocity  in  the  physical  model.  Displacement  velocity  for  the 
downbound  tows  was  far  less  consistent.  Some  of  the  tows 
created  an  upstream  velocity  ( Steve  Kuhr  test  06)  or  a  reduction 
in  the  ambient  velocity  ( Steve  Kuhr  test  08).  Both  of  these 
results  were  anticipated,  based  on  observation  of  the  physical 
model.  But  in  some  of  the  downbound  tests,  the  displacement 
velocities  were  in  the  downbound  direction  ( John  Matthews  tests 
04  and  10) ,  which  is  opposite  to  observations  in  the  physical 
model.  This  may  have  been  the  result  of  the  bow  effect  from  the 
towboat.  In  the  Harold  Turner  three-wide  downbound  tests 
(tests  9,  11,  and  13),  the  displacement  velocities  approached 
zero,  which  was  radically  different  from  the  upbound  tests  and 
from  observations  in  the  physical  model. 

c.  Displacement  velocity  near  stern  of  barge(s).  This  location  was 
selected  to  define  the  combined  effect  of  displacement  velocity, 
suction  effects  of  the  propellers,  and  the  bow  wave  at  the  front 
of  the  towboat  for  tests  with  one  unloaded  barge. 

d.  Velocities  entering  the  propellers.  These  velocities  are  the 
result  of  the  suction  effects  in  the  propellers,  the  displace¬ 
ment  effects  of  the  towboat,  and  the  wake  effects  of  the  barges. 
Only  suction  and  displacement  effects  occur  for  the  tests  having 
one  barge  ahead  of  the  towboat,  and  both  effects  are  in  a  direc¬ 
tion  opposite  the  tow  direction  based  on  observation  of  the 
physical  model.  In  the  upbound  tests,  the  velocities  are  down¬ 
stream  for  most  of  the  tests  as  expected.  In  the  downbound 
tests,  many  of  the  tests  have  velocities  in  an  upbound  direc¬ 
tion.  But  in  several  of  the  downbound  tests  ( Steve  Kuhr  test 
08;  John  Matthews  tests  02  and  04;  Harold  Turner  tests  02,  04, 
06,  08,  11,  13),  velocities  are  in  the  same  direction  as  the 
towboat,  which  is  opposite  to  the  expected  direction  or  the 
direction  observed  in  the  upbound  tests.  It  should  be  noted 
that  the  velocities  beneath  the  towboat  were  the  maximum  ob¬ 
served  at  any  of  the  five  locations  in  the  prototype  tests.  It 
was  also  noted  that  dye  injection  in  the  flow  visualization 
tests  showed  strong  vertical  components  near  the  propellers,  and 
two-dimensional  velocity  measurements  were  suspect  in  this 
region. 

e.  Velocities  in  propeller  wash.  These  velocities  are  the  result 
of  two  opposing  velocities.  The  propeller  jet  acts  in  a  direc¬ 
tion  opposite  to  the  tow  direction.  The  wake  flow  velocities 
act  in  the  same  direction  as  the  tow.  Most  of  the  downbound 
tests  had  bottom  velocities  in  a  downstream  direction,  which  is 
opposite  that  observed  in  the  physical  model.  Only  John 
Matthews  tests  02  and  04  had  propeller  wash  velocities  opposite 
the  tow  direction.  For  the  upbound  tows,  about  half  of  the  tows 
had  downstream  propeller  wash,  but  magnitudes  were  less  than 
0.3  fps.  For  the  other  upbound  tows,  the  propeller  wash  was 
skewed  in  a  generally  upstream  direction. 
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48.  Based  on  data  from  this  set  of  prototype  tests,  the  variations  be¬ 
tween  flow  direction  in  model  and  prototype  suggest  three  possibilities: 

a.  The  flow  direction  indicated  by  the  prototype  meters  is  incor¬ 
rect. 

b.  The  understanding  of  velocities  around  moving  tows  is  incorrect. 

c.  The  tow  was  not  centered  over  the  meter  and  the  wake  flow  was 
measured,  not  the  propeller  v-ash. 
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PART  VII:  DEVELOPMENT  OF  PREDICTIVE  RELATIONS 


49.  Using  model  and  prototype  data,  predictive  relations  were  developed 
for  the  following  three  areas: 

a.  Return  velocity  for  vessels  sailing  off  the  channel  center  linf 
in  nonprismatic  channels . 

b.  Propeller  jet  velocities  for  moving  tows. 

c.  Velocities  immediately  beside  and  beneath  the  vessel  due  to 
displacement  effects . 

In  these  relations  the  following  sign  convention  will  be  used:  (a)  downstream 
flows  are  positive,  (b)  upstream  flows  are  negative,  and  (c)  for  slack— water 
conditions,  flows  opposite  tow  direction  are  positive.  All  velocities  are  as¬ 
sumed  to  act  parallel  to  the  axis  of  the  channel,  which  is  the  same  as  the 
tow,  and  tow-induced  and  ambient  velocities  are  added  or  subtracted,  depending 
on  direction. 


Return  Velocities 


Method  1:  Effective  Area  Method 

50.  Two  methods  of  computing  return  velocity  were  developed  in  this 
study . 

51.  Return  velocities  were  determined  in  this  analysis  using  the 
Bouwmeester  et  al.  (1977)  method.  Modifications  were  required  to  determine 
the  effective  area  for  vessels  sailing  at  any  distance  off  the  channel  center 
line  and/or  in  a  natural  channel  cross  section.  Marchal  and  Spronck  (1977) 
treated  each  side  of  the  vessel  as  being  independent  of  the  other  for  the  ves¬ 
sel  sailing  off  the  channel  center  line.  Using  this  assumption,  the  waterway 
area  used  in  the  Bouwmeester  equation  would  be  twice  the  area  of  the  side  for 
which  return  velocity  is  being  computed.  This  assumption  was  tested  using  the 
Olmsted  model  data  shown  in  Plates  10-12.  The  average  return  velocities  were 
computed  using  the  Bouwmeester  method  for  the  left  and  right  sides  of  the 
vessel  using  effective  area  2  x  Aside  and  effective  width  2  x  BBide  and  are 
shown  in  the  following  tabulation: 
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Tow 

Location 

Computed 

V 

r ,  Left 

f  DS 

„  2A 

N  —  Left 

Left  A 

A 

m 

Computed 

V 

r ,  Right 

f  DS 

M  2A 

N  =  Right 

Right 

A 

m 

A 

0.35 

45 

0.40 

37 

B 

0.31 

51 

0.45 

31 

C 

0.19 

69 

1.1 

13 

Note:  ALe£t  -  waterway  area  on  left  side  of  vessel. 

A,,,  =  submerged  cross-sectional  area  of  midship  section. 


The  comparison  with  the  model  data  is  good  with  the  exception  of  Vr  Riaht  at 
tow  location  C.  The  measured  return  velocities  are  less  than  the  computed, 
indicating  that  a  portion  of  the  return  flow  on  the  right  side  is  either  going 
under  the  vessel  or  passing  to  the  left  (in  front)  of  the  vessel.  One  way  to 
handle  this  is  to  define  an  effective  area  Aeff  and  effective  width  Beff 
that  should  be  used  when  the  vessel  sails  close  to  a  bank  line.  PIANC  (1987) 
presents  a  method  for  determining  an  effective  area  for  vessels  sailing  off 
the  center  line,  but  this  method  is  applicable  only  to  prismatic  channels.  A 
method  is  needed  that  can  define  an  effective  area  and  width  for  any  shape  of 
channel  when  the  vessel  sails  near  the  bank.  The  empirical  relation  used 
herein  is 

Aeff  *  A^to,  (2ABide)  (49> 

where 


Ae££  —  effective  area  for  determining  return  velocity  for  vessels 
sailing  near  a  bank  line 

Aside  =  area  of  side  for  which  return  velocity  is  being  computed 
Afacf.or  =  [A0/(2Aai,ie)  ]  P 

P  -  empirical  coefficient 
The  effective  width  is  determined  by 


Beff  =  Bfactot  (2Bside) 


(50) 


where 


B#f£  -  effective  width  for  determining  return  velocity  for  vessels 
sailing  near  a  bank  line 
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^factor  **  [B0/(2B8ide)]p 


Bside  -  surface  width  of  side  for  which  return  velocities  are  being 
computed 

Note  that  the  Afactor  and  Bfaccor  equal  1  for  a  vessel  sailing  on  the  center 
line  of  a  prismatic  channel.  The  average  depth  is  computed  from 


h  = 


>f  £ 


Je£f 


(51) 


Blockage  ratio  and  then  return  velocity  are  computed  for  each  side  of  the  ves¬ 
sel.  To  determine  P  .  it  is  known  that  P  -  0  for  all  Naide  >  31  ,  where 
Nside  is  the  blockage  ratio  for  each  side  of  the  vessel,  based  on  the  results 
in  the  tabulation.  For  VR  right  at  location  C  of  the  Olmsted  data,  various 
P  values  were  tested  and  P  =  0.475  resulted  in  good  agreement  between  ob¬ 
served  and  calculated  values.  Plotting  these  data  in  Plate  16  with  P  versus 
Nsid,,  defines  a  tentative  relationship  for  determination  of  Aef{  and  Beff 
for  vessels  sailing  in  channels  having  low  NBide  .  More  data  are  needed  to 
refine  this  relationship,  but  good  correlation  was  found  using  this  effective 
area  to  compute  Schijf's  limiting  velocity  and  the  method  developed  by  Fuehrer 
and  Romisch  (1977)  given  by  Equation  44. 

Method  2 :  total  area  method 

52.  One  of  the  problems  associated  with  the  effective  area  method  is 
that  in  some  cases  the  return  velocity  equations  show  that  the  ship  is  travel¬ 
ing  at  speeds  greater  than  the  critical  speed  when  the  tow  is  close  to  one 
bank.  When  this  happens,  no  solution  is  obtained.  In  the  total  area  method, 
the  average  return  velocity  is  computed  for  the  entire  cross-section  and  then 
this  value  is  proportioned  on  the  port  and  starboard  sides  of  the  vessel 
(Vca)  depending  on  the  position  of  the  tow  in  the  cross  section.  This  tech¬ 
nique  eliminates  almost  all  of  the  critical  speed  problems  associated  with 
method  1.  Another  problem  with  method  1  involves  the  use  of  the  Bouwmeester 
equation  in  an  asymmetric  channels.  One  of  the  required  inputs  for  Bouw¬ 
meester  is  the  bank  slope,  which  can  be  difficult  to  define  in  an  asymmetric 
channel.  For  that  reason,  method  2  used  the  original  Schijf  (1949)  equation 
(without  correction)  since  it  requires  only  h  ,  N  ,  and  V  .  Using  data 
shown  in  Table  8  from  the  Delit  Hydraulics  Laboratory  (1979)  and  from 
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1:33-  and  i: /O-scale  models  reported  in  Maynord*  the  following  equation  was 
obtained 


Zif  =  0.42  skew  +  0.58  (52) 


where  skew  =  AQ/(2Aside) .  Method  2  was  used  to  compute  the  Olmsted  model  con¬ 
ditions  given  in  Plates  10—12  as  follows: 


Tow 

Computed 

V 

r 

Location 

Left 

Right 

A 

0.35 

0.38 

B 

0.32 

0.40 

C 

0.26 

0.61 

Agreement  with  the  observed  data  is  good  but  additional  data  are  needed. 

53.  Once  the  average  return  velocity  for  each  side  of  the  vessel  is 
determined,  the  distribution  from  vessel  to  bank  is  needed.  Using  data  from 
the  Delft  Hydraulics  Laboratory  (1979)  and  the  1:35-  and  1: 70-scale  models 
(Table  8),  the  ratio  of  maximum  return  velocity  near  the  vessel  Vrsm  to  the 
average  return  velocity  for  the  side  of  the  vessel  was  found  to  have  the 
following  relationship: 


a 


Iff  =  0 . 024  (Nside)  -  0.75 


(53) 


where  N  =  side 

side 

A 

m 

This  relationship  is  similar  to  Fuehrer  and  Romisch  (1977)  except  that  an  area 
ratio  (Nside)  is  used  rather  than  a  width  ratio.  The  ratio  a  was  found  to  be 
a  good  parameter  for  defining  the  type  of  velocity  distribution  as  follows: 


*  S.  T.  Maynord.  1989  (May).  "Return  Velocity  Distribution  and  Flow  Visual¬ 
ization  for  Commercial  Navigation,"  letter  report,  US  Army  Engineer  Water¬ 
ways  Experiment  Station,  Vicksburg,  MS. 
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q  Value 

Distribution 

Equation 

a  <  1 . 0 

Uniform 

vr(y)  - 

Q  >  1.0 

Exponential 

vt(y) 

=  a  exp 

-c 

(  v  Y1 

y  -  b 

vr8 

"B.ide  -  b 

J  , 

where  C  =  2(a  -  1)  .  Additional  data  are  needed  to  better  define  the 
relationship  of  C  and  a  . 


Propeller  Jet  Velocities 


54.  Both  the  Verhey  (1983)  and  Oebius  (1984)  methods  for  propeller  jet 
velocities  were  evaluated  using  some  of  the  physical  model  results.  The 
Verhey  method,  which  includes  the  wake  effects  found  to  be  significant  in  this 
study,  gave  the  most  realistic  results.  For  these  reasons  the  Verhey  method 
was  chosen  for  further  testing  and  development  using  the  physical  model  data. 
The  following  presents  several  modifications  and  how  input  parameters  were 
defined  in  the  Verhey  method: 

a.  The  following  propeller  thrust  coefficients  at  zero  speed  of 
advance  from  Prosser  (1986)  are  used  in  the  Verhey  method: 

Pitch/propeller  diameter 

Kto  (Kort  nozzle) 


0.8  0.9  1.0 
0.37  0.44  0.51 


b.  The  Verhey  method  was  calibrated  to  velocities  2  ft  above  bottom 
only . 

c.  The  wake  coefficient  was  defined  as 


Maximum  of 


W  =  0.4 


*  0.04 


or 


W  =  0. 3 


(54) 


where  Dep  is  the  local  depth  of  flow.  Draft  and  beam  are  for 
either  the  barges  or  the  towboat,  whichever  has  the  maximum 
cross-sectional  area.  This  results  in  values  of  W  ranging 
from  0.3  to  0.52  for  typical  drafts  <  9  ft,  depth  >  15  ft,  and 
beam  <  105  ft . 
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d.  Advance  coefficient  Ja  is  specified  <  1  and  is  used  to  modify 
the  thrust  coefficient  at  zero  speed  of  advance  according  to 


Kt  =  Kto(l  -  J.) 


(55) 


e.  The  Verhey  method  is  modified  to  account  for  two  propellers  and 
uses  distance  from  each  propeller  to  determine  velocity  from 
each  propeller  and  then  adds  the  two  velocities  together  as 
shown  in  Plate  17.  The  combined  velocities  are  limited  to  the 
maximum  velocity  in  the  jet. 

f.  Verhey  uses  the  contraction  diameter  D0  to  define  the  strength 
of  the  wake  flow  at  the  channel  bottom.  A  better  representative 
size  would  depend  on  the  vessel  cross  section  and  the  water 
depth.  Conceptually  the  strength  of  the  wake  at  the  bottom 
behind  the  vessel  will  increase  for  increasing  vessel  cross 
section  and  will  decrease  for  increasing  clearance  beneath  the 
vessel.  Several  relations  were  tried  and  the  following  empiri¬ 
cal  relation  was  found  to  give  satisfactory  characteristic 
dimension  of  the  wake  flow  strength: 


DW  =  0 . 1 


.  bd  1 

[Dep  -  d J 


(56) 


Beam  and  draft  are  for  either  the  barges  or  towboat,  whichever 
has  the  maximum  cross-sectional  area. 

g.  The  coefficients  C  and  b'  are  used  in  the  Verhey  method  to 
compute  the  propeller  jet  contribution  to  the  total  velocity. 

The  coefficient  C  varies  the  location  where  the  maximum  jet 
velocity  occurs,  and  the  coefficient  b  varies  the  magnitude 
and  rate  of  decay  downstream  of  the  propellers.  Verhey  used  the 
relation  (VQ  -  V)/VQ  to  reduce  C  for  the  moving  tow. 

Results  from  the  physical  model  demonstrated  a  better  correla¬ 
tion  with  a  parameter  frequently  used  in  propeller  jet  studies, 
distance  from  propeller  axis  to  channel  bottom/propeller  diam¬ 
eter.  The  empirical  relation  developed  for  C  is 


C  =  0. 17  -  0. 39  ^ 
Dp 
or 

C  =  0.04 


4  3 


Maximum  of 


(57) 


This  C  was  used  in  Equations  21-23.  The  coefficient  b'  was 
found  to  vary  with  advance  coefficient  Ja  and  h^/D^ 
according  to 


b' 


0.52  +  J,  +  0.24 


hp 


(58) 


Variation  of  C  and  b'  with  hp/Dp  accounts  for  the  limited 
jet  spreading  that  occurs  because  of  the  channel  bottom. 


55.  A  comparison  of  velocities  calculated  with  the  Verhey  method  with 
physical  model  velocities  is  shown  in  Plates  18-25.  The  Verhey  method  was  not 
compared  to  the  prototype  data  due  to  the  uncertainty  of  the  prototype  tow 
location . 

Displacement  Velocity  Near  Vessel 

56.  Velocities  generated  at  the  bow  of  the  vessel  Vfcb  in  the  same 
direction  as  the  vessel  were  used  in  the  analysis  of  the  following  equation 


Vbb  =  /(V,  Dep,  d,  b) 


When  grouped  into  dimensionless  parameters 


(59) 


* »  -v 

^bb  f  d  b 

V  Dep  ’  Dep 


(60) 


Channel  size  and  return  velocity  effects  were  found  to  be  small  when  data  from 
the  three  physical  model  channel  sizes  were  compared  and  were  not  used  in  the 
analysis.  The  resulting  best-fit  equation  based  on  the  physical  model  data 
for  bow  velocity  is 


0.18 


'10.6 


Dep 


Dep 


1.6* 


(61) 
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Observed  and  predicted  prototype  and  model  bow  velocities  are  listed  in 
Tables  9  and  10,  respectively.  Because  of  the  bow  effects  of  the  towboat  for 
the  tests  with  one  unloaded  barge,  the  only  prototype  tests  that  can  be  com¬ 
pared  to  the  model  derived  regression  equations  are  Harold  Turner  tests  05—13. 
The  observed  prototype  values  (Table  9)  compare  reasonably  well  with  computed 
values  from  the  model  derived  equation. 

57.  Even  in  extremely  wide  waterways  having  essentially  zero  average 
return  velocity,  tow  movement  can  create  significant  displacement  velocities 
immediately  beside  and  beneath  the  barges  and  towboat.  The  magnitude  of  these 
velocities  will  be  dependent  primarily  on  the  water  depth,  tow  draft,  beam, 
and  vessel  speed.  The  analysis  of  Fuehrer  and  Romisch  (1977)  in  Equation  45 
was  tested,  but  results  showed  that  the  beam  width  is  also  an  important  param¬ 
eter,  even  in  wide  waterways.  A  comparison  of  bottom  velocities  from  similar 
vessels  in  the  unconfined,  one  near— bank,  and  confined  channels  is  shown  in 
Table  11.  The  velocity  immediately  beneath  the  vessel  Vbd  at  the  axis  of 
the  vessel  is  relatively  independent  of  channel  size  ana  return  velocity.  The 
Fuehrer  and  Romisch  equation  also  shows  this  due  to  the  limited  influence  of 
1/N  in  Equation  45.  A  regression  analysis  of  the  physical  model  data  was 
conducted  using  the  following  dimensionless  ratios 


V  J 


Dep  ’  Dep 


(62) 


and  resulted  in  the  following  equation 


0.16 


(63) 


Observed  and  predicted  prototype  and  model  velocities  are  shown  in  Tables  12 
and  13,  respectively.  Similar  to  the  bow  velocities,  the  prototype  displace¬ 
ment  velocities  observed  for  the  one  unloaded  tow  cannot  be  compared  to  the 
model  relations  because  of  the  towboat  influence.  Again  only  the  Harold 
Turner  tests  05-13  can  be  compared  to  the  model-derived  relations  for  dis¬ 
placement  velocity.  In  every  case,  the  model  relations  overestimated  the 
center-line  displacement  velocities  observed  in  the  prototype.  These 
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differences  between  model  and  prototype  velocities  are  likely  caused  by  a 
combination  of  the  following  factors: 

a.  The  velocity  meter  was  15  in.  above  bottom  in  the  prototype  and 
24  in.  above  the  bottom  in  the  model. 

b.  The  effects  of  ambient  flows  are  unknown.  All  model  tests  were 
conducted  in  slack  water.  Prototype  upbound  and  downbound  tests 
exhibited  significant  differences.  In  this  analysis  the  ambient 
flows  were  added  or  subtracted  from  the  observed  velocity 
depending  on  tow  direction. 

c.  The  model  velocity  meter  had  a  threshold  velocity  of  0.5  fps , 
which  means  that  model  regression  analysis  did  not  have  data 
below  this  value  while  much  of  the  prototype  data  were  below 
this  value. 

d.  The  bottom  roughness  was  greater  in  the  prototype.  This  will 
move  the  higher  velocities  away  from  the  channel  bottom. 

e.  Most  importantly,  the  prototype  meter  did  not  have  a  very  fast 
frequency  response  and  was  not  capturing  the  peak  velocity  that 
occurs  near  the  bow  of  the  barges.  For  this  reason  the  model- 
derived  relation  is  retained  for  estimating  velocities  in  the 
prototype . 

58.  Velocities  immediately  adjacent  to  the  moving  vessel  will  also  be 
significant  even  in  waterways  having  negligible  average  return  velocity.  Com¬ 
parison  of  the  velocities  measured  27.5  ft  from  the  edge  of  the  barges  (Vout. 
in  Table  11)  shows  that  the  channel  size  is  a  significant  parameter,  or  stated 
otherwise,  that  the  return  velocity  is  contributing  to  the  total  Vout  .  Vout 
is  similar  to  Vrmax  used  by  Fuehrer  and  Romisch  (1977)  with  one  important 
difference.  Vout  is  a  bottom  velocity  whereas  Vr  max  is  a  depth-averaged 
velocity.  It  is  likely  that  the  deptn-averaged  velocity  and  the  bottom  veloc¬ 
ity  at  2  ft  are  relatively  close  because  the  boundary  layer  has  not  grown. 
Determination  of  Vout  will  require  the  following  type  of  analysis: 


Vout  =  /(displacement  velocity) 


r.  /return  velocity  or' 
3  \blockage  ratio 


(64) 


which  is  similar  to  the  approach  given  by  Fuehrer  and  Romisch  (1977)  for  ve¬ 
locities  beneath  the  vessel.  The  problem  with  this  approach  is  that  model  and 
prototype  data  are  not  available  to  evaluate  this  approach  fully  Initially, 
this  approach  was  discarded  and  the  following  relation  used  by  Fuehrer  and 
Romisch  (1977)  was  evaluated: 
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(65) 


The  correlation  between  these  two  parameters  was  very  poor.  An  approximation 
was  used  that  allowed  evaluation  of  Equation  64.  The  same  factors  driving 
Vbd  in  Equation  63  were  assumed  to  drive  the  /(displacement  velocity)  terr.i 
in  Equation  64.  Data  were  needed  that  had  a  negligible  contribution  from  the 
term  /(blockage  ratio)  in  Equation  64.  The  only  data  meeting  this  require¬ 
ment  were  the  Steve  Kuhr  and  John  Matthews  prototype  data  that  were  discussed 
in  Part  VI.  From  this  prototype  data  the  coefficient  in  the  /(displacement 
velocity)  term  was  derived 
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Dep 


0.54 


/(blockage  ratio) 


(66) 


Note  that  (b/Dep)°  5<‘(d/Dep)0'68  was  taken  from  the  equation  for  Vbd  .  The 
next  step  in  the  analysis  was  to  subtract  the  first  term  of  the  right  side  of 
Equation  66  from  the  observed  Vout/V  for  all  of  the  physical  model  data  and 
the  Harold  Turner  prototype  data.  This  left  the  term  /(blockage  ratio)  , 
which  was  evaluated  as  a  function  of  the  blockage  ratio  N  .  The  resulting 
equation  is 
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Observed  and  predicted  model  and  prototype  velocities  are  compared  in 
Table  14.  Computed  prototype  values  were  generally  larger  than  observed  and 
may  have  been  a  result  of  uncertainty  in  the  location  of  the  tow  relative  to 
the  meter  as  well  as  other  factors  mentioned  in  paragraph  57. 

Relationship  Between  Velocity  and  Shear  Stress 

59.  The  forr e/unit  area  or  shear  stress  is  often  a  more  descriptive 
parameter  than  bottom  or  depth-averaged  velocity.  For  a  given  shear  stress, 
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bottom  velocity  varies  widely  depending  on  distance  from  the  bottom.  Most 
sediment  transport  relations  are  based  on  shear  stress.  Only  a  few  are  based 
on  velocity,  and  most  of  these  use  depth— averaged  velocity.  Shear  stress  has 
rarely  been  easily  or  accurately  estimated  from  a  single  point  velocity  or 
depth-averaged  velocity.  Relations  (to  be  discussed  subsequently)  exist  to  do 
this,  but  the  problem  lies  in  their  application.  With  most  open-channel  flow 
problems,  the  boundary  layer  is  fully  developed  and  extends  from  the  bottom  to 
the  water  surface.  For  vessel-induced  flows,  the  boundary  layer  has  not  be¬ 
come  fully  developed,  which  makes  the  determination  of  shear  stress  more  dif¬ 
ficult  and  the  use  of  velocity  less  reliable.  The  following  paragraphs  pre¬ 
sent  ^^th^dc  for  determining  the  relationship  between  shear  stress  and  the 
velocities  induced  by  the  vessel: 

a.  Return  velocities  from  vessel  to  bank.  This  method  (Blaauw 
et  al.  1984)  uses  the  equation 


r 


(68) 


where 

r  -  shear  stress 

C£r  =  local  skin  friction  coefficient 
p  =  water  density 

The  local  skin  friction  coefficient  is  defined  by  Schlicting 
(1968)  as 
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ft 


(2.87 


1.58  log 


5 


(69) 


where 

X  =  distance  from  beginning  of  boundary  layer 
development  to  maximum  velocity 

ks  *»  equivalent  sand  roughness 

=  3  to  4  d50 

Blaauw  et  al.  (1984)  define  X  as 


X  = 


vt  -  v 


L 


(70) 
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where  L  is  the  ship  length.  Equations  68-70  provide  a  method 
of  transferring  return  velocity  into  shear  stress  acting  on  the 
channel  bottom. 

b.  Propeller  iet  velocity.  In  the  same  method  as  in  subpara¬ 
graph  a,  Blaauw  and  van  de  Kaa  (1978)  define  Cfr  =  0.06  -  0.11 
for  propeller  jet  flow.  In  a  later  study  Verhey  (1983)  says 
that  Cfr  can  vary  from  0.06  to  0.36.  The  higher  values  of  Cfr 
in  the  propeller  jet  versus  the  return  flow  are  due  to  the 
higher  turbulence  levels.  Bottom  velocity  from  the  propeller 
jet  equations  is  used  in  Equation  68.  The  distance  from  the 
beginning  of  the  propeller  jet  velocity  to  the  maximum  propeller 
jet  velocity  is  approximately  four  times  the  clearance  from  the 
bottom  of  the  propeller  to  the  channel  bottom.  A  local  skin 
friction  coefficient  of  two  to  three  times  the  value  obtained 
from  Equation  69  is  recommended  for  propeller  jet  flows. 

c.  Displacement  velocities  beneath  vessel.  The  distance  from  the 
flow  reversal  at  the  bow  of  the  lead  barge  to  the  point  of  maxi¬ 
mum  velocity  due  to  the  displacement  current  is  three  to  four 
times  the  clearance  beneath  the  bottom  of  the  barge.  This  dis¬ 
tance  is  used  for  X  ,  and  Equations  68  and  69  are  proposed  for 
determining  the  shear  stress  caused  by  the  displacement  flow 
beneath  the  vessel. 


49 


PART  VIII:  SUMMARY  AND  CONCLUSIONS 


60.  Analysis  of  the  existing  literature  demonstrated  thr  i.->1 1  -  _ng  needs 
relative  to  navigation  effects  in  large  waterways  found  in  many  parts  of  the 
United  States: 

a.  Propeller  jet  (wash)  velocity  is  strongly  affec<~"d  by  tow  speed, 
and  any  predictive  relation  should  cons  id*- ■  speed  explicitly. 

The  width  (if  any)  of  propeller  jet  attack  on  the  channel  bottom 
needs  to  be  defined. 

b.  Return  velocity  relations  need  to  be  applicable  to  vessels  mov¬ 
ing  in  asymmetric  channels  at  various  positions  in  the  cross 
section  with  various  ambient  velocities. 

c.  The  use  of  velocity  in  navigation  effects  studies  needs  to  be 
standardized,  and  methods  are  needed  for  converting  velocity  to 
bed  shear  stress. 

61.  Using  flow  visualization  and  model  and  prototype  measurements,  this 
investigation  has  demonstrated  the  following  flow  patterns  near  the  path  of  a 
moving  tow: 

a.  Bow  velocities.  Bow  velocities  are  dependent  primarily  on  the 
water  depth,  draft,  beam  width,  and  vessel  speed.  Bow  veloc¬ 
ities  act  in  the  same  direction  as  the  vessel  and  have  been 
recorded  as  far  as  550  ft  in  front  of  the  vessel  for  three- 
barge-wide  loaded  tows  in  21  ft  of  water.  Maximum  bow  bottom 
velocities  occurred  just  beneath  the  bow  of  the  vessel  and  were 
up  to  2.9  fps  for  a  three-barge-wide  loaded  tow  traveling  in 

15  ft  of  water. 

b.  Center-line  displacement  velocites.  Center-line  displacement 
velocity  Vbd  under  the  barges  or  the  towboat  is  also  primarily 
dependent  on  water  depth,  draft,  beam  width,  and  vessel  speed. 

Vbd  acts  opposite  the  direction  of  the  tow  and  reached  a  maxi¬ 
mum  at  a  location  three  to  four  times  the  clearance  beneath  the 
vessel  downstream  from  the  bow.  Maximum  Vbd  exceeded  4  fps 
for  a  three-barge-wide  loaded  tow  traveling  in  15  ft  of  water. 

c.  Velocities  beneath  towboat.  Velocities  in  this  region  are 
affected  by  many  factors  including  wake  flow  from  barges,  flow 
entering  propellers,  and  displacement  velocity  from  the  towboat. 
Flow  near  the  propellers  is  highly  three-dimensional,  and 
neither  the  model  or  prototype  velocity  measurements  can  be 
considered  reliable  in  this  zone.  Flow  visualization  demon¬ 
strated  that  the  yarn  strings  and  dye  used  as  indicators  pointed 
straight  up  as  the  propellers  passed  over  the  indicator 
position . 

d.  Propeller  iet  velocities.  Velocities  in  the  propeller  wash 
region  are  affected  by  many  factors  including  wake  flow  from 
barges  and  towboat,  propeller  type,  size  and  speed,  open  wheel 
or  Kort  nozzle,  water  depth,  and  tow  speed.  For  the  moving  tow, 
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these  factors  combine  to  produce  a  relatively  narrow  propeller 
jet  path  (60  ft  wide)  moving  opposite  to  the  tow  and,  on  each 
side  of  the  propeller  jet,  a  wake  flow  moving  in  the  same  direc¬ 
tion  as  the  tow.  The  60-ft  width  of  propeller  jet  needs  better 
definition  regarding  the  effects  of  such  factors  as  tow  speed, 
towboat  size,  and  propeller  speed  and  size. 

e.  Velocities  outside  path  of  vessel.  Unlike  the  velocities  be¬ 
neath  the  vessel,  velocity  outside  the  path  of  the  vessel  was 
significantly  affected  by  channel  size,  meaning  that  the  return 
velocity  was  contributing  to  the  total  velocity. 


62.  Predictive  relations  were  developed  for  return  velocities,  propeller 
jet  velocities,  and  displacement  velocities.  For  vessels  sailing  off  the 
channel  center  line  and/or  in  nonprismatic  channels,  two  methods  were  devel¬ 
oped  so  that  average  return  velocity  could  be  computed  for  each  side  of  the 
vessel.  A  method  was  proposed  for  estimating  the  distribution  of  return  cur¬ 
rents  from  vessel  to  shore.  Additional  studies  are  needed  to  define  the  in¬ 
fluence  of  ambient  currents.  Photographic  techniques  were  demonstrated  as 
being  an  effective  tool  for  determining  the  distribution  of  return  velocity. 

63.  The  Verhey  (1983)  method  for  propeller  jet  velocities  was  applied 
using  the  physical  model  data.  Several  empirical  coefficients  required  in  the 
Verhey  method  were  determined  from  the  physical  model  data.  The  primary  ad¬ 
vantages  of  this  method  are  the  incorporation  of  vessel  speed  and  the  wake 
flow  effects  in  determining  propeller  jet  velocity. 

64.  After  the  appropriate  physical  parameters  were  selected,  empirical 
relations  were  developed  for  bow  velocity,  displacement  velocity  beneath  the 
vessel,  and  velocity  outside  the  path  of  the  vessel. 

65.  Relations  were  presented  for  determining  the  relation  between  veloc¬ 
ity  and  shear  stress  that  are  based  on  accepted  boundary  layer  theory. 
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Table  1 

Bottom  Velocity  Tests 
Unconfined  Channel 


Velocity 


Test 

Propeller 

Barge 

Draft 

Depth 

Speed 

Meter 

No. 

Sneed,  mm 

Configuration* 

ft 

ft 

fos 

Location** 

15 

185 

3W 

x  2L 

3 

15 

11.1 

\ 

14 

135 

3W 

*  2L 

15 

6.8 

\ 

22 

185 

3W 

x  2L 

30 

11.1 

\ 

21 

135 

3W 

x  2L 

30 

6.7 

\ 

13 

185 

1W 

x  1L 

15 

11.2 

B 

12 

135 

1W 

x  1L 

15 

6.6 

B 

29 

185 

1W  x  1L 

30 

11.5 

B 

28 

135 

1W 

x  1L 

30 

6.7 

B 

9 

185 

3W  x  2L 

2 

15 

11.1 

\ 

8 

135 

3W  x  2L 

15 

6.5 

\ 

25 

185 

3W  x  2L 

30 

11.4 

\ 

24 

135 

3W 

x  2L 

30 

6.8 

\ 

11 

185 

1W 

x  1L 

15 

11.0 

B 

10 

135 

1W  x  1L 

15 

6.7 

B 

27 

185 

1W 

x  1L 

30 

11.6 

B 

26 

135 

1W 

x  1L 

30 

6.8 

B 

16 

TOW 

3W 

x  2L 

3 

15 

10.6 

l 

\ 

23 

TOW 

3W 

x  2L 

30 

11.1 

18 

185 

3W 

x  1L 

15 

11.3 

17 

135 

3W 

x  1L 

15 

6.8 

20 

185 

3W 

x  1L 

30 

11.3 

19 

135 

3W 

x  1L 

30 

6.5 

31 

185 

3W  x  2L 

30 

11.5 

z 

30 

135 

30 

6.8 

2 

39 

185 

15 

11.1 

- 

38 

135 

15 

6.5 

2 

40 

TOW 

15 

11 . 1 

2 

42 

1 8  5  f 

15 

8.1 

\ 

41 

1 3  5  f 

15 

5.9 

\ 

Notes:  All  tests  conducted  with  Kort  nozzle. 

Propeller  rotation:  When  viewed  from  the  stern,  left  propeller  has 
counterclockwise  rotation;  right  propeller  has  clockwise  rotation. 

*  W  -  width:  L  =  length. 

**  A  =  center  line  and  80  ft  from  center  line,  both  2  ft  above  bottom. 
B  ==  center  line  and  45  ft  from  center  line,  both  2  ft  above  bottom 
C  =  26  and  45  ft  from  center  line,  both  2  ft  above  bottom, 
f  No  towing  assist. 


Table  2 

Bottom  Velocity  Tests 
Near-Bank  Channel 


Velocity 

Test 

Propeller 

Barge 

Draft 

Depth 

Speed 

Meter 

No . 

Speed,  rom 

Configuration* 

ft 

ft 

f  D  S 

Location** 

201 

185 

3W  x  2L 

( 

9 

15 

11.2 

t 

\ 

200 

135 

3W  x  2L 

15 

6.7 

i 

209 

185 

3W  > 

<  2L 

30 

11.4 

t 

208 

135 

3W  x  2L 

30 

6.8 

i 

\ 

203 

185 

1W  x  1L 

15 

11.3 

B 

202 

135 

1W  x  1L 

15 

6.7 

B 

211 

185 

1W  ) 

<  1L 

30 

11.3 

B 

210 

135 

1W  ; 

<  1L 

30 

6.8 

B 

228 

185 

3W  > 

<  2L 

2 

15 

11.2 

l 

\ 

227 

135 

3W  ) 

<  2L 

15 

6.5 

l 

220 

185 

3W  > 

<  2L 

30 

11.4 

l 

\ 

219 

135 

3W  3 

<  2L 

30 

6.7 

t 

\ 

207 

185 

1W  x  1L 

15 

11.4 

B 

206 

135 

1W  x  1L 

15 

6.8 

B 

213 

185 

1W  ) 

<  1L 

30 

11.2 

B 

212 

135 

1W  ) 

<  1L 

30 

6.5 

B 

205| 

185 

3W  ) 

<  2L 

< 

9 

15 

11.2 

A 

204f 

135 

15 

6.7 

i 

3 

218| 

185 

30 

11.2 

t 

\ 

2 1 7 1 

135 

30 

6 . 6 

l 

\ 

229 

185 

15 

11.2 

( 

-* 

230 

135 

15 

6.6 

231 

TOW 

15 

12.2 

215 

185 

30 

11.2 

214 

135 

30 

6 . 7 

216 

TOW 

30 

6 . 7 

114 

TOW 

15 

10.5 

A 

117 

TOW 

30 

11.4 

l 

\ 

222  ft 

150 

3W  ; 

<  1L 

21 

12  . 1 

A 

221ft 

100 

3W  ; 

<  1L 

8 . 5 

A 

226ft 

150 

1W  ) 

<  1L 

18.1 

B 

225ft 

100 

1W  : 

<  1L 

8.8 

B 

224ft 

150 

1W  x  1L 

2 

18.1 

B 

223ft 

100 

1W  ; 

<  1L 

l 

13.3 

B 

Notes : 

Propeller 

rotation:  When  viewed 

from 

the  stern, 

left  prope 

ller 

has 

clockwise  rotation;  right  propeller  has  counterclockwise  rotation. 


*  W  =  width;  L  =  length. 

**  A  =  center  line  and  80  ft  from  center  line,  both  2  ft  above  bottom. 
B  =  center  line  and  45  ft  from  center  line,  both  2  ft  above  bottom. 
C  =  26  and  45  ft  from  center  line,  both  2  ft  above  bottom. 

|  Open-wheel  test. 

It  Reproduces  prototype  tests  of  MV  Harold  Turner. 


Table  3 

Bottom  Velocity  Tests 
Confined  Channel 


Test 

No . 

Propeller 
Sneed,  rpm 

Barge 

Confi duration* 

Draft 

ft 

Depth 

ft 

Speed 

f  DS 

Velocity 

Meter 

Location** 

301 

185 

3W  : 

x  2L 

( 

9 

15 

10.0 

l 

\ 

300 

135 

3W  : 

x  2L 

15 

6.6 

A 

321 

185 

3W  x  2L 

30 

11.2 

i 

\ 

320 

135 

3W  : 

x  2L 

30 

6.6 

i 

\ 

310 

185 

1W  : 

x  1L 

15 

11.3 

B 

309 

135 

1W  x  1L 

15 

6.6 

B 

333 

185 

1W  : 

x  1L 

30 

11.0 

B 

332 

135 

1W  : 

x  1L 

30 

6.3 

B 

304 

185 

3W  : 

x  2L 

2 

15 

10.6 

t 

\ 

303 

135 

3W  x  2L 

15 

6.6 

i 

\ 

327 

185 

3W 

x  2L 

30 

11.2 

i 

\ 

326 

135 

3W 

x  2L 

30 

6.6 

l 

\ 

307 

185 

1W 

x  1L 

15 

11.3 

B 

306 

135 

1W 

x  1L 

15 

6.5 

B 

330 

185 

1W 

x  1L 

30 

11.3 

B 

329 

135 

1W 

x  1L 

30 

6.5 

B 

3 1 6  f 

185 

3W 

x  2L 

9 

15 

11.5 

L 

\ 

3 1 5 1 

135 

15 

6.6 

L 

\ 

3 1 8  f 

185 

30 

11.2 

A 

3 1 7  f 

135 

30 

6.6 

l 

\ 

313 

185 

15 

11.2 

( 

-J 

312 

135 

15 

6.7 

314 

TOW 

15 

12.0 

324 

185 

30 

11.2 

323 

135 

30 

6.6 

325 

TOW 

30 

12.9 

302 

TOW 

15 

9.7 

t 

\ 

319 

TOW 

30 

12.6 

l 

\ 

Notes:  Propeller  rotation:  When  viewed  from  the  stern,  left  propeller  has 

clockwise  rotation;  right  propeller  has  counterclockwise  rotation. 

*  W  =  width;  L  =  length. 

**  A  -  center  line  and  80  ft  from  center  line,  both  2  ft  above  bottom. 
B  -  center  line  and  45  ft  from  center  line,  both  2  ft  above  bottom. 
C  -  26  and  45  ft  from  center  line,  both  2  ft  above  bottom, 
f  Open-wheel  test. 


Table  4 

Steve  Kuhr  Prototype  Tests 


Data  File  for  Towboat  Study:  summer  1987 
Date:  29  July  1987 

Towboat  name:  Steve  Kuhr  HP  Rating:  1,500 

Towboat  dimensions  (L  x  W  x  D) ,  ft:  70  x  24  x  7 
Number  of  screws:  3  Open  versus  Kort:  Open 

Barge  Configuration:  1  wide  x  1  long  (Unloaded) 

Barge  Dimension  (L  x  W  x  D) :  195  x  35  x  1.25 

Overall  Tow  Dimensions  (L  x  W) :  265  x  35 

Propeller  Diameter,  in.:  60  Propeller  Pitch,  in.:  46 

Red  Gear  Ratio:  5.17:1 

Approximate  River  Velocity,  fps :  0.46 


Run 

Up  (U) 
or 

Down  CD) 

Propeller 
Soeed.  mm 

Speed,  fps 

Clock  ( 

L1072901 

U 

155 

6.6 

1028 

02 

D 

155 

7.4 

1051 

03 

U 

230 

9.5 

1110 

04 

D 

230 

11.0 

1124 

05 

U 

310 

13.9 

1137 

06 

D 

310 

14.7 

1149 

07 

U 

155 

7.4 

1439 

08 

D 

155 

7.8 

1500 

09 

U 

230 

9.8 

1515 

10 

D 

230 

10.6 

1533 

11 

U 

310 

13.9 

1546 

12 

D 

310 

13.9 

1556 

13 

U 

310 

14.7 

1610 

Results 

Runs  01  to  06  in  shallow  water,  about  16-17  ft  deep.  Vessel  exhibited  notice¬ 
able  squat  (about  9  in.)  only  at  highest  speed/rpm. 

Runs  07  to  12  in  deeper  water,  about  21-22  ft  deep.  Vessel  exhibited  minimal 
squat  at  all  speeds/rpms. 


Speed  indicated  is  relative  to  land  tracking  station. 


Table  5 

John  Matthews  Prototype  Tests 


Data  File  for  Towboat  Study:  summer  1987 
Date:  31  July  1987 

Towboat  name:  John  Matthews  HP  Rating:  4,200 

Towboat  dimensions  (L  x  W  x  D) ,  ft:  148  x  35  x  8.5 
Number  of  screws:  2  Open  versus  Kort:  Kort 

Barge  Configuration:  1  wide  X  1  long  (Unloaded) 

Barge  Dimension  (L  x  W  x  D) :  195  x  35  x  1.25 

Overall  Tow  Dimensions  (L  x  W) :  335  x  35 

Propeller  Diameter,  in.:  96  Propeller  Pitch,  in.:  84 

Red  Gear  Ratio:  3.958:1 

Approximate  River  Velocity,  fps :  0.29 


Up  (U) 


Down  (D) 


Propeller 
Speedy  rpn 


Speed,  fps 


Clock  (EDT) 


L1073101 

U 

100 

8.4 

1107 

02 

D 

100 

8.6 

1123 

03 

U 

160 

13.3 

1134 

04 

D 

160 

13.8 

1145 

05 

U 

225 

15.7 

1203 

06 

D 

225 

20.3  (18.8)* 

1220 

07 

U 

100 

8.0 

1531 

08 

D 

100 

8.4 

1550 

09 

U 

160 

13.3 

1605 

10 

D 

160 

13.3 

1617 

11 

U 

225 

19.2 

1633 

12 

D 

225 

19.2 

1648 

*  (18.8)  based  on  event  marks. 

Results 

Runs  01  to  06  in  shallow  water,  about  16-17  ft  deep.  Vessel  exhibited  notice¬ 
able  squat  (up  to  about  24  in.)  only  at  highest  speed/rpm. 

Runs  07  to  12  in  deeper  water,  about  21-22  ft  deep.  Vessel  exhibited 
considerable  squat  (up  to  about  18  in.)  at  highest  speeds/rpm. 

Speed  indicated  is  relative  to  land  tracking  station. 


Table  6 

Harold  Turner  Prototype  Tests 


Data  File  for  Towboat  Study:  summer  1987 


Date:  5  August  1987 

Towboat  name:  Harold  Turner  HP  Rating:  5,600 


Towboat  dimensions  (L  x  W  x  D) ,  ft:  140  x  42  x  9 


Number  of  screws :  2 

Propeller  Diameter,  in.:  110 
Red  Gear  Ratio:  4.345:1 
Approximate  River  Velocity,  fps : 


Open  versus  Kort:  Kort 

Propeller  Pitch,  in.:  100 

0.50 


Run 

Up  (U) 
or 

Down  (D) 

Propeller 

Speed, 

rom 

Speed , 

f  D  S 

Clock 

(EDT'I 

Loaded  (L) 
or 

Unloaded  fUl 

Overall 
Dimensions , 
ft 

L1080501 

U 

100 

13.4 

1120 

1W  x 

1L  U 

335L  x 

42W 

02 

D 

100 

12.9 

1130 

1W  x 

1L  U 

335L  x 

42W 

03 

U 

150 

16.7 

1154 

1W  x 

1L  U 

335L  x 

42W 

04 

D 

150 

16.7 

1208 

1W  x 

1L  U 

335L  x 

42W 

05 

U 

100 

9.6 

1323 

1W  X 

1L  L 

335L  x 

42W 

06 

D 

100 

9.3 

1339 

1W  x 

1L  L 

335L  x 

42W 

07 

U 

150 

15.9 

1400 

1W  x 

1L  L 

335L  x 

42W 

08 

D 

150 

17.6 

1413 

1W  x 

1L  L 

335L  x 

42W 

09 

D 

165 

8.3 

1520 

3W  x 

5L  L 

1 , 115L  x 

105W 

10 

U 

100 

9.0 

1545 

3W  x 

1L  L 

335L  x 

105W 

11 

D 

100 

9.6 

1607 

3W  x 

1L  L 

335L  x 

10  5W 

12 

U 

150 

12.4 

1627 

3W  x 

1L  L 

335L  x 

105W 

13 

D 

150 

12.9 

1645 

3W  x 

1L  L 

335L  x 

105W 

Results 

Runs  01  to  13  all  in  deep  water,  about  22  ft  deep,  and  farther  from  Kentucky 
bank  than  any  other  runs.  Run  09  is  F.  M.  Baker,  an  opportunistic  tow  that 
had  an  identical  towboat  to  Harold  Turner. 


Speed  indicated  is  relative  to  land  tracking  station. 


Water  Ambient  _ Observed  Velocity,  fps** _ 

Depth  Tow  Upbound  (U)  Propeller  Bottom  _ Displacement _ 

at  Tow  Config-  Test  or  Speed  V  Velocity  Bow  of  Bow  of  Stern  of  Entering  Propeller 

ft  urat ion*  No .  Downbound  (D)  rpro  fps  fps  Barge  Barge  Barge  Propeller  Wash 
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Table  8 

Return  Velocity  Data 


Identifier 

V 

ri 

V 

r 

Side 

Skew 

N 

side 

V 

r  em 

V 

r  s 

Delft  Hydraulics 

Laboratory  (1979) 

Figure  6 

0.83 

Port 

0.57 

63.5 

2.34 

2.18 

Starboard 

3.96 

9.2 

1.00 

Figure  7 

0.79 

F'>rt 

0.57 

63.5 

2.26 

2.49 

Starboard 

3.96 

9.2 

1.00 

Figure  5 

1.00 

Port 

1.00 

36.3 

1.71 

1.00 

Starboard 

1.00 

36.3 

1.71 

1:35  scale* 

Test  22 

1.32 

Port 

1.50 

12.6 

1.06 

0.84 

Starboard 

0.75 

25.2 

1.29 

Test  23 

1.60 

Port 

3.00 

6.3 

1.00 

0.87 

Starboard 

0.60 

31.5 

1.47 

Test  24 

0.94 

Port 

1.00 

18.9 

1.13 

1.05 

Starboard 

1.00 

18.9 

1.05 

1:70  scale* 

Test  47,  47A 

1.00 

Port 

1.00 

12.5 

1.05 

1.00 

Starboard 

1.00 

12.6 

1.07 

Test  48,  48A 

0.94 

Port 

0.66 

19.0 

1.03 

1.17 

Starboard 

2.03 

6.2 

1.00 

Test  54,  54A 

1.10 

Port 

1.43 

8.8 

1.00 

0.95 

Starboard 

0.7/ 

16.3 

1.15 

*  S.  T.  Maynord.  1989  (May).  "Return  Velocity  Distribution  and  Flow 

Visualization  for  Commercial  Navigation,"  letter  report,  US  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  MS. 


Table  9 

Observed  Versus  Computed  Prototype  Bow  Velocities 


Test 

No. 

Tow 

Config¬ 

uration* 

Upstream 

or 

Downstream 

(U) 

(D) 

Propeller 

Speed 

rom 

V 

f  DS 

Ambient 

Bottom 

Velocity 

f  DS 

V 

bb 

Observed 

,  fps* 

Compared 

05 

11L 

U 

100 

9.6 

0. 

,  15 

0.19(119) 

0.47(360) 

06 

11L 

D 

100 

9.3 

0.76(168) 

0.69(180) 

07 

11L 

U 

150 

15.9 

?? 

08 

11L 

D 

150 

17.6 

1.0(142) 

1.2(180) 

09 

35L 

D 

165 

8.3 

1.4(144) 

1.1(180) 

10 

31L 

U 

100 

9.0 

1.2(94) 

0.96(360) 

11 

31L 

D 

100 

9.6 

1.3(148) 

1.2(180) 

12 

31L 

U 

150 

12.4 

1.1(90) 

1.4(360) 

13 

31L 

D 

150 

12.9 

2.0(119) 

1.6(180) 

Note:  All  data  are  from  Harold  Turner  tow. 

*  11U  =  ]  wide  x  1  long  unloaded 

11L  ■  1  wide  x  1  long  loaded 

35L  =  3  wide  x  5  long  loaded 

31L  -  3  wide  x  1  long  loaded 

**  Value  in  parentheses  following  velocity  is  azimuth,  with  a  value  of 
180  deg  being  downstream. 


15 

2 

.90 

2 

.77 

207 

0 

.00 

0 

.12 

14 

1 

.90 

1 

.70 

206 

0 

.00 

0 

.07 

22 

0 

.50 

0 

.58 

205 

3 

.20 

2 

.80 

21 

0 

.00 

0 

.  35 

204 

1 

.  80 

1 

.67 

13 

1 

.  55 

1 

.44 

218 

0 

.  57 

0 

.  59 

12 

1 

.00 

0 

.85 

217 

0 

.00 

0 

.  35 

29 

0 

.00 

0 

.31 

114 

2 

.30 

2 

.62 

11 

0 

.00 

0 

.12 

117 

0 

.50 

0 

.60 

10 

0 

.00 

0 

.07 

222 

2 

.  75 

1 

.42 

27 

0 

.00 

0 

.02 

221 

1 . 

.23 

1 

.00 

16 

2 

.50 

2 

.65 

226 

1. 

.60 

1 

.  10 

23 

0 

.00 

0 

.  58 

225 

0. 

.00 

0 

.  53 

18 

4 

.  10 

2 

.82 

224 

0, 

.00 

0. 

.09 

17 

1 

.55 

1. 

.70 

223 

0. 

.00 

0. 

.06 

20 

0 

.57 

0 

.59 

301 

2. 

.  15 

2, 

.  50 

19 

0 

.00 

0. 

.34 

300 

1. 

.60 

1. 

,65 

42 

1 

.55 

2. 

.02 

321 

0. 

00 

0. 

.59 

41 

1. 

.40 

1. 

.47 

320 

0. 

00 

0. 

35 

201 

3 

.45 

2. 

.80 

310 

1 . 

.30 

1. 

46 

200 

1. 

80 

1. 

67 

309 

0. 

57 

0. 

85 

209 

0. 

57 

0. 

60 

333 

0. 

00 

0. 

30 

208 

0. 

.00 

0. 

36 

304 

0. 

00 

0. 

22 

203 

1. 

40 

1. 

46 

315 

1. 

55 

1. 

65 

202 

0. 

80 

0. 

86 

318 

0. 

00 

0. 

59 

211 

0. 

00 

0. 

30 

317 

0. 

00 

0. 

35 

210 

0. 

00 

0. 

18 

302 

1. 

55 

2. 

42 

220 

0. 

00 

0. 

05 

319 

0. 

00 

0. 

66 

Table  11 

Comparison  of  Model  Bottom  Velocities  for  Three  Channel  Types  (9-ft  Draftl 


Ves 

;sel 

Spe 

?ed . 

fos 

Te 

st  No  . 

V  CL 

bd 

• 

fps 

V 

out 

.  ff 

is 

Width 

Depth 

Uncon 

INear 

Conf 

Uncon 

!N°«r 

Conf 

Uncon 

INear 

Conf 

Uncon 

INear 

Conf 

3W 

15 

11. 

.  1 

11. 

,2 

10. 

0 

15 

201 

301 

3 

.  7 

4. 

1 

3 

.  9 

2 

.  2 

2 

.6 

3 

.4 

15 

6. 

.8 

6. 

.7 

6. 

6 

14 

200 

300 

2. 

.  7 

2. 

9 

2, 

.4 

1 

.3 

1 

.  5 

2 

.0 

30 

11. 

.  1 

11. 

.4 

11. 

2 

22 

209 

321 

1. 

.2 

1. 

5 

1 

.7 

1 

.  3 

1 

.4 

1 

.6 

30 

6 

.  7 

6 

.8 

6. 

6 

21 

208 

320 

0 

.7 

0. 

8 

0 

.8 

0 

.0 

0 

.0 

1 

.  5 

1W 

15 

11 

.2 

11 

.3 

11. 

3 

13 

203 

310 

2 

.6 

2. 

7 

2 

.4 

1 

.2 

1 

.  2 

1 

.2 

15 

6. 

.6 

6. 

.  7 

6. 

6 

12 

202 

309 

1 

.4 

1. 

2 

1 

.2 

0 

.6 

0 

.0 

1 

.4 

30 

11. 

5 

11. 

3 

11. 

0 

29 

211 

333 

0. 

.  7 

0. 

7 

0. 

.  7 

0 

.0 

0 

.0 

1 

.  1 

3W 

15 

10 

,6 

10 

.  5 

9. 

7 

16* 

114 

302 

2. 

.3 

3. 

7 

3 

.0 

1 

,3 

2 

.3 

2 

.5 

30 

11 

.1 

11. 

A 

12. 

6 

23* 

117 

319 

1 

.4 

1. 

5 

1 

.  7 

1 

.2 

1 

.4 

1 

.  5 

V  26 

V 

45 

— 

bd 

— 

— 

bd 

— 

3W 

15 

11. 

.  1 

11. 

.2 

11. 

2 

39 

229 

313 

3 

,2 

3. 

5 

3 

.4 

2 

.4 

2. 

.8 

2 

4 

15 

6. 

.  5 

6 

.6 

6. 

7 

38 

230 

312 

2 

.0 

2. 

9 

3 

.2 

1 

.7 

2 

.  3 

2 

.6 

30 

11. 

.  5 

11 

.2 

11. 

2 

31 

215 

324 

1 

.8 

1. 

5 

1 

.  5 

1 

.6 

1 

.  5 

1 

.2 

30 

6. 

.8 

6 

.  7 

6. 

6 

30 

214 

323 

1 

.2 

0. 

8 

1 

.  1 

1 

,0 

0 

,6 

1 

.  1 

15 

11 

.  1 

12 

.2 

12. 

0 

40* 

231 

314 

2 

,  8 

3. 

2 

3 

.2 

2 

.2 

2 

.2 

2 

.4 

Note:  Uncon  -  unconfined 

INear  -  near-bank 

Conf  -  confined 

Vb<JCL  -  bottom  displacement  velocity  at  tow  centei  line. 

Vbd26  -  bottom  displacement  velocity  at  26  ft  from  tow  center  line. 

Vbd45  -  bottom  displacement  velocity  at  45  ft  from  tow  center  line. 

Vout  “  bottom  displacement  velocity  at  27.5  ft  from  edge  of  barge. 
*  Towing  test. 


Table  12 

Observed  Versus  Computed  Prototype  Center-Line  Displacement  Velocities 


Test 

No . 

Tow 

Config¬ 

uration* 

Upstream 

or 

Downstream 

(U) 

(D) 

Propeller 

Speed 

rnm 

V 

fps 

Ambient 

Bottom 

Velocity 

f  D  S 

V 

bb 

Observed 

,  fps* 

Compared 

05 

11L 

U 

100 

9.6 

0. 

,15 

0.71(209) 

1.3(180) 

06 

11L 

D 

100 

9.3 

0. 14(?) 

0.95(360) 

07 

11L 

U 

150 

15.9 

0.8(200) 

2.0(180) 

08 

11L 

D 

150 

17.6 

0 . 47 ( ?) 

1.9(360) 

09 

35L 

D 

165 

8.3 

0.09(101) 

1.6(360) 

10 

31L 

U 

100 

9.0 

1.4(120) 

2.1(180) 

11 

31L 

D 

100 

9.6 

0.09(135) 

1.9(360) 

12 

31L 

U 

150 

12.4 

1.4(129) 

2.8(180) 

13 

31L 

D 

150 

12.9 

0.12(88) 

2.6(360) 

Note:  All  data  are  from  Harold  Turner  tow. 

*  1.1U  =  1  wide  x  1  long  unloaded 

11L  «*  1  wide  x  1  long  loaded 

35L  =  3  wide  x  5  long  loaded 

31L  -=  3  wide  x  1  long  loaded 

**  Value  in  parentheses  following  velocity  is  azimuth,  with  a  value  of 
180  deg  being  downstream. 


Table  13 

Observed  Versus  Computed  Model  Center-Line  Displacement  Velocities 


Test 

V 

,  fps 

Test 

V 

,  fps 

tb 

bb 

No. 

Observed 

Computed  I 

No. 

Observed 

Computed 

15 

3. 

.  70 

3. 

.58 

207 

1. 

,15 

0. 

.73 

14 

2. 

,70 

2. 

.19 

206 

0. 

65 

0. 

,43 

22 

1. 

,20 

1. 

.54 

205 

4. 

,20 

3. 

.62 

21 

0. 

,  70 

0. 

.  93 

204 

2. 

80 

2. 

.  16 

13 

2. 

.60 

1, 

.99 

218 

1. 

.80 

1, 

.  55 

12 

1. 

,40 

1. 

.17 

217 

1. 

10 

0. 

,  91 

29 

0. 

,70 

0, 

.87 

114 

3. 

.70 

3. 

.39 

11 

0. 

,60 

0 

.70 

117 

1, 

.50 

1. 

.  58 

10 

0. 

,50 

0. 

.42 

222 

4. 

.00 

2. 

.59 

27 

0. 

,50 

0 

.32 

221 

2. 

.20 

1. 

.82 

16 

2. 

,  30 

3. 

.42 

226 

2. 

,30 

2. 

.14 

23 

1. 

.40 

1 

.54 

225 

0, 

.90 

1. 

.04 

18 

5. 

.80 

3. 

.65 

224 

0. 

.65 

0. 

.77 

17 

2. 

.70 

2 

.19 

223 

0, 

.65 

0. 

.56 

20 

1 

.70 

1 

.56 

301 

3. 

.90 

3. 

.23 

19 

0 

.80 

0 

.90 

300 

2. 

.40 

2, 

.13 

42 

3, 

.  10 

2 

.61 

321 

1. 

.70 

1, 

.55 

41 

1 

.55 

1 

.90 

320 

0. 

.75 

0. 

.91 

201 

4. 

.10 

3 

.62 

310 

2. 

.40 

2. 

.01 

200 

2. 

.90 

2. 

.16 

309 

1. 

20 

1. 

.17 

209 

1. 

,50 

1. 

.58 

333 

0. 

74 

0. 

.84 

208 

0. 

80 

0. 

.94 

304 

1. 

00 

1. 

23 

203 

2. 

.70 

2. 

.01 

315 

0. 

78 

2. 

13 

202 

1. 

.15 

1. 

.19 

318 

1. 

.55 

1. 

.55 

211 

0. 

,70 

0. 

.86 

317 

1. 

06 

0. 

91 

210 

0. 

.50 

0. 

,52 

302 

3. 

00 

3. 

13 

220 

0. 

.50 

0. 

.  57 

319 

1. 

70 

1. 

74 

Table  14 


Observed  Versus  Computed  V  ,  Model  and  Prototype 
_ out _ 


Physical  Model  Data _  _ Prototype  Data 


Test 

No. 

N 

V 

out 

Test 

No. 

N 

V 

out 

Observed 

Computed 

Observed 

Computed 

15 

15. 

.00 

2. 

,22 

2. 

.26 

1 

67. 

9 

0.30 

0. 

42 

14 

15, 

.00 

1 , 

.32 

1 

.  38 

2 

67  . 

9 

0.15 

0. 

40 

22 

27. 

.60 

1. 

.35 

1, 

,07 

3 

67. 

9 

0.30 

0. 

59 

21 

27  , 

.60 

0. 

.00 

0 

.65 

4 

67. 

9 

0.16 

0. 

62 

13 

38. 

.  10 

1 . 

23 

1. 

.  11 

5 

67. 

9 

0.71 

0. 

85 

12 

38, 

.10 

0. 

.62 

0, 

.65 

6 

67. 

9 

0.52 

0. 

84 

29 

76. 

.20 

0. 

.00 

0 

.51 

7 

134. 

5 

0.22 

0. 

30 

11 

171. 

.40 

0. 

.97 

0 

.34 

8 

134. 

5 

0.17 

0. 

28 

10 

171. 

.40 

0. 

,00 

0 

.20 

10 

134. 

5 

0.22 

0. 

39 

27 

342. 

.90 

0. 

.00 

0 

.16 

12 

134. 

5 

0.29 

0. 

52 

16 

15. 

.00 

1. 

.26 

2, 

.16 

13 

134. 

5 

0.37 

0. 

59 

23 

27. 

.60 

1. 

.23 

1 

.07 

1 

45. 

0 

0.08 

0. 

73 

18 

15. 

.00 

2. 

,  50 

2 

.30 

3 

45. 

0 

0.78 

1. 

15 

17 

15. 

.00 

1. 

.75 

1. 

.38 

4 

45. 

0 

0.75 

1. 

14 

20 

27 

.60 

0. 

.85 

1. 

.09 

5 

45. 

0 

1.77 

1. 

30 

19 

27. 

.60 

0. 

.00 

0, 

.63 

7 

67. 

9 

0.50 

0. 

49 

42 

15 

.00 

1. 

.47 

1. 

.65 

8 

67. 

9 

0.30 

0. 

48 

41 

15. 

.00 

0. 

.00 

1 

.20 

9 

67. 

9 

0.50 

0. 

81 

201 

11 

.10 

2. 

.57 

2 

.60 

10 

67. 

9 

0.40 

0. 

77 

200 

11 

.10 

1 

.49 

1 

.56 

11 

67. 

9 

0.69 

1. 

16 

209 

21 

.80 

1. 

.43 

1 

.24 

12 

67. 

9 

0.84 

1. 

13 

208 

21 

.80 

0 

.00 

0 

.74 

1 

112. 

7 

0.73 

0. 

81 

203 

25 

.00 

1 

.25 

1 

.31 

1 

112. 

7 

0.59 

0. 

72 

202 

25 

.00 

0. 

.00 

0 

.78 

3 

112. 

7 

0.99 

1. 

00 

211 

38 

.  10 

0 

.00 

0 

.  69 

4 

112. 

7 

0.78 

0. 

95 

210 

38 

.10 

0 

.00 

0 

.41 

5 

112. 

7 

0.55 

0. 

59 

220 

57 

.10 

0 

.00 

0 

.  46 

6 

112. 

7 

0.86 

0. 

51 

207 

85 

.  70 

1 

.00 

0 

.43 

7 

112. 

7 

1 . 19 

0. 

96 

206 

85 

.  70 

0 

.00 

0 

.26 

8 

45. 

1 

1.53 

1. 

80 

205 

11 

.10 

2 

.83 

2 

.60 

9 

45. 

1 

1.22 

0. 

82 

204 

11 

.  10 

1. 

.47 

1 

.  56 

10 

45. 

1 

0.90 

1. 

OC 

218 

21 

.80 

0 

.98 

1 

.22 

11 

45. 

1 

0.93 

0. 

96 

217 

21 

.80 

1 

.26 

0 

.  71 

12 

45. 

1 

1.44 

1. 

36 

114 

11 

.10 

2 

.29 

2 

.44 

13 

45. 

1 

1.89 

1 . 

31 

117 

21 

.80 

1 

.43 

1 

.  24 

— 

— 

— 

- 

— 

222 

15 

.60 

2 

,28 

1 

.93 

— 

— 

— 

- 

— 

221 

15 

.60 

1 

.61 

1 

.36 

— 

— 

— 

- 

— 

226 

29 
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Photo  3.  Surface  current  patterns,  1 : 70-scale  tow,  location  C 


Photo  4.  Surface  current  patterns,  1: 120-scale  tow,  location 


Surface  current  patterns,  1: 120-scale  tow,  location 


Photo  6.  Surface  current  patterns,  1: 120-scale  tow,  location 
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NOTE:  HORIZONTAL  AXIS  IS  DISTANCE.  TOW  LENGTH  DEFINES  THE  SCALE 
FOR  THIS  AXIS.  TOW  LENGTH  FOR  2  BARGES  PLUS  TOWBOAT  IS  599  FT 
TOW  LENGTH  FOR  1  BARGE  PLUS  TOWBOAT  IS  404  FT. 


SMALLEST  DIVISION  FOR  VELOCITY  METER  READING  IS  5  UNITS.  LARGE 
DIVISION  FOR  VELOCITY  METER  READING  IS  25  UNITS.  RELATIONSHIP 
BETWEEN  VELOCITY  METER  READING  AND  PROTOTYPE  VELOCITY  IS  GIVEN  IN 
i ABlE  A  1 
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1'20-SCALE  MODEL 
TESTS  22  AND  23 
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1  Ml  7715 


TEST  39 


BOTTOM  VELOCmES 

1'80-SCALE  MODEL 
TESTS  31  AND  39 
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BOTTOM  VELOCITIES 

1'20-SCALE  MDDEL 
TESTS  207  AND  208 
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PLATE  A1 3 


TEST  223 


BOTTOM  VELOCmES 

1‘20-SCALE  MDDEL 
TESTS  222  AND  223 
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1  .'20-SCALE  MODEL 
TESTS  224  AND  225 
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TEST  318 
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TESTS  317  AND  318 
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BOTTOM  VELOCmES 

1.20-SCALE  MODEL 
TESTS  319  AND  320 
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TEST  329 
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1'20-SCALE  MdDEL 
TESTS  321  AND  329 
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a'  Eccentricity  factor  for  vessel  sailing  off  channel  center  line 
a, A  Coefficients 

A  Function  of  propeller  height  above  bottom  and  if  rudder  is 
present  behind  propeller 

a,B  Coefficients  in  Hochstein  (1967)*  return  velocity  equations 

Ae££  Effective  waterway  area  for  determining  return  velocity  for 
vessel  sailing  near  bank 

A,„  Submerged  cross-sectional  area  of  midship  section 
Ac  Cross-section  area  of  waterway 

Aaid«  Waterway  area  on  one  side  of  vessel  for  which  return  velocity  is 
being  computed 

b'  Coefficient 

b  Beam  of  vessel 

Bc  Surface  width  of  waterway 

Be££  Effective  waterway  width  for  determining  return  velocity  for 
vessel  sailing  near  bank 

Baide  Surface  width  of  side  for  which  return  velocities  are  being 
computed 

C  Coefficient 

C'  Coefficient 

C£t  Local  skin  friction  coefficient 
d  Draft  of  vessel 

D0  Orifice  diameter  or  contraction  diameter 
Dp  Propeller  diameter 

D'  Modified  propeller  diameter  (Oebius  1984) 

Dep  Local  depth  of  flow 

DW  Characteristic  dimension  of  wake  flow 

DX  Distance  from  propeller  center  to  maximum  velocity  (Oebius 
1984) 


*  References  cited  in  this  appendix  are  included  in  the  References  at  the  end 
of  the  main  text. 
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E  Coefficient 


F(o),K2 
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hp 
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J. 
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K 

Kt 

Kto 
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N 

Np 


NBid» 
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Coefficients  in  return  velocity  distribution  equation 
Acceleration  due  to  gravity 
Average  channel  depth  =  A0/B0 

Distance  from  center  line  of  propeller  to  bottom 
Modified  advance  coefficient 
Advance  coefficient 
Equivalent  sand  roughness 
Constrainment  factor  for  determining  Vcr 
Thrust  coefficient  for  moving  navigation 
Thrust  coefficient  at  zero  speed  of  advance 
Vessel  length 

Propeller  speed,  revolutions  per  second 
Blockage  ratio  -  A0/A,„ 

Number  of  propellers 

Blockage  ratio  for  each  side  of  the  vessel 
Coefficient  for  determing  A„ff  and  Beff  ;  engine  power 


J2  -f-  2 

Z  y 

R£  Blade  radius  (from  outside  of  hub  to  blade  tip) 
R„  Hub  radius 

s  Cotangent  of  side  slope  angle 
U0  Ambient  velocity  in  undisturbed  channel 
U2  Velocity  increment 
V  Vessel  speed  relative  to  earth 
Vs  Entrance  or  advance  velocity 


Vb(x)  Bottom  velocity  in  x-direction  at  coordinate  x 


Vbb 

Vbd 


Vb 


,max 
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b, 


max, J=0 
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cr 


ve 

Vi 

Vo 


v„ 
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vr(y) 

v„ 
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ram 


vw 

vx 

V(r) 


W 

x 


x 


O 


X 

y 

z 

zl 

Z 


Bottom  velocity  at  bow  of  vessel  acting  in  the  same  direction 
Displacement  velocity  beneath  vessel 

Maximum  bottom  velocity  for  moving  navigation  (v  /  0) 

Maximum  bottom  velocity  for  maneuvering  navigation  (V  =  0) 
Limiting  velocity  of  self-propelled  vessel 
Bottom  velocity  in  the  propeller  wash  region 

Velocity  at  infinity  assumed  equal  to  vessel  speed  (Oebius  1984) 

Orifice  velocity  at  outlet  or  jet  velocity  at  propeller 

Bottom  displacement  velocity  27.5  ft  outside  the  edge  of  the 
barge 

Return  velocity 

Maximum  return  velocity  at  the  vessel 

Variation  of  return  velocity  with  distance  from  vessel 
Average  return  velocity  for  each  side  of  the  vessel 
Maximum  return  velocity  near  tow 
Wake  velocity 

Velocity  in  x-direction  at  coordinates  x,r 

Maximum  propeller  jet  velocity  at  radial  distance  r  from 
propeller  axis 

Velocity  in  x-direction  at  coordinates  x,r  -  o 
Wake  fraction 

Distance  from  outlet  or  propeller  measured  along  jet  axis 
Limit  of  flow  establishment  zone 

Distance  from  beginning  of  boundary  layer  development  to 
maximum  velocity 

Horizontal  distance  from  the  center  line  of  propeller 

Water-level  drawdown 

Distance  from  canal  axis  to  vessel 

Vertical  coordinate  measured  from  center  line  of  propeller 
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a  Ratio  of  maximum  to  average  return  velocity 
aQ  Fuehrer  and  Romisch  (1977)  parameter  for  flow  beneath  vessel 
Coefficient  (Oebius  1984) 
p  Water  density 

t  Shear  stress 
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